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FOREWORD 


The technology program described herein was sponsored 
by the National Aeronautics and Space Administration, 
Lyndon B. Johnson Space Center, Houston, Texas, under 
Contract NAS9-12077. The study was conducted during 
the 15-month period from 29 June 1972 to 28 September 
1973. The NASA technical monitor was Mr, R. C. Kahl. 

At Rocketdyne, Mr. L. P. Combs was program manager and 
Dr. C. L, Oberg was project engineer. 

The literature review, Priem-type stability analysis, 
and steady-state combustion model analysis were done 
primarily by Dr. R. C. Kesselring. The NREC model 
analysis was done by Mr. M. D. Schuman. The generalized 
cavity analysis was done by Dr. C. Warner III. Overall 
technical direction was provided by Dr. C. L. Oberg. 


R-94S3 

ii 



ABSTRACT 


Analytical studies have been made of the relative com- 
bustion stability of various propellant combinations 
of interest to NASA-JSC when used with hardware con- 
figurations representative of current design practices 
and with or without acoustic cavities. Available li- 
terature was reviewed to locate and summarize available 
experimental results relating to stability. Two com- 
bustion instability models, a Priera-type model and a 
modification of the Northern Research and Engineering 
(NREC) instability model, were used to predict the var- 
iation in engine stability with changes in operating 
conditions, hardware characteristics or propellant com- 
bination, exclusive of acoustic cavity effects. The 
NREC model was developed for turbojet engines but is 
applicable to liquid propellant engines. A steady-state 
combustion model was used to predict the needed input 
for the instability models. In addition, preliminary 
development was completed on a new model to predict the 
influence of an acoustic cavity with specific allowance 
for the effects the nozzle, steady flow and combustion. 
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NOMENCLATURE 


ARABIC SYMBOLS 

A = specific acoustic admittance of nozzle 

A^ = cross-sectional area of chamber 

Ap = Priera stability index (critical disturbance amplitude) 

A^ = cross-sectional area of nozzle throat 

= iterative expansion coefficient, defined by Eq. 66 
Cp = orifice discharge coefficient 

C^ = droplet drag coefficient 

Cp * . thrust coefficient 

Cp = heat capacity — constant pressure 

Cy = heat capacity--constant volume 

CR = contraction ratio 

= NREC combustion coefficients (i = 1, 2z, 2r, 29, 3, 4, 5, 6) 
c = isentfopic sound velocity 

c* = characteristic exhaust velocity 

D = chamber diameter; with subscript, droplet diameter 

Dp = mass median diameter of fuel spray 

D ) = mass median diameter of oxidizer spray 

ox ^ ' 

~ mass diffusivity 

E = energy content used in NREC model 

e^ = total internal energy (thermal only) of spray 

e^ = energy release rate (see Eq. 3^ 

F = force of interaction between spray and gas 

F = thrust 
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S(r|T^) 


J ( ) 


1 . 2,3 


thrust per element 
Green's function 

modified Green's function defined by Eq. 68 

2 

gravitational constant (32.2 Ibm-ft/lbf-sec ) 
heat of combustion 
heat of vaporization 

source coefficient for generalized cavity damping model, Eq. 57 

Bessel function of first kind and order m 

parameters used in Priem-type model, Eq, 4, 5, and 6 
B/c 

thermal conductivity of gas 

chamber length 

depth of acoustic cavity 

Priem burning rate parameter, Eq. 1 

Mach number 

defined by Eq. 8 

mass accumulation parameter, Eq. 3 
mixture ratio 
molecular weight 

derivative of burning rate, Ae/Az 
unit normal vector directed out of volume 
Nusselt number for heat transfer 
Nusselt number for mass transfer 
oscillatory source term, Eq, i6 
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X 



Pr 

P 

P 

Pc 

Pv 

Q 

Qs 

R 

Re 


s 

S(?;p) 

T 

t 

t 

w 

u 

u 

u 

s 

V 


= Prandtl number, uC /k 

P g 

= pressure 

= oscillatory pressure 
= steady-state chamber pressure 
= vapor pressure 

= heat release per unit volume from homogeneous reactions 
= heat release per unit volume of chamber from spray combustion 
= gas constant; also, chamber radius 
= Reynold’s number 

= response factor from Nusselt number, Eq. 48 
= response factor,^ "^tiu 

- response factor from heat blockage term, Eq. 53 
* radial position 
= chamber radius 
= position vector 
= area 

= source terra defined by Eq. 64 
= temperature 
= flame temperature 

= expansion coefficient defined by Eq, 80 
= time 

= Priem wave time, 2irR/c • 

= velocity 

= oscillatory velocity 
= spray velocity 
= volume 
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gas 

V, 

drops 

AV 

V. . 
inj 

w(r) 

w 


w. 


inj 


w 


vap 


OX 


y 

z 


axial gas velocity 
axial spray velocity 

Priem parameter, (''gas - ''drops>^= 
liquid injection velocity 
weight factor, Eq. 77 

mass generation of gas per unit chamber volume from spray combustion 
injected mass flowrate 

mass generation of gas per unit chamber volume from spray combustion 
burned fraction of fuel 
burned fraction of oxidizer 

specific acoustic admittance, defined by Eq. 59 
axial coordinate 


GREEK SYMBOLS 


a 

6 

Y 

6(r-r') 
6 . . 

6 ( ) 
e 


0 

A 

A 

P 


damping coefficient, imarginary part of Q 
complex angular frequency, 
heat capacity ratio, Cp/C^ 

Dirac delta function 

Kronecker delta (5^^ = 1 ifi = j, = 0 ifi^j) 

variational operator 

overall burned fraction of spray 

eigenvalue for closed, rigid walled chamber, Eq. 62 

angular coordinate 

normalization factor, Eq. 63 

parameter in source distribution, A = ,u T 

gas density 
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Pit 

a 

X 


OJ 


w. . 

inj 


= spray density, mass per unit chamber volume 
= density-diameter product for droplet 
= fu-u j w 
= NREC delay time 

= complex eigenvalue, <t> - kr^ = wr /c+jar^/c; also with 
subscript, eigenfunction 

= angular frequency 

= Priem parameter, w. . = w. . m/A 
^ ’ inj inj c 


SUPERSCRIPTS 


— (overbar) 


(i) 

A 

SUBSCRIPTS 

0 
c 
d 

dr 

f 

g 

1 

N 

m 

q 

ox 

r 

s 

t 

z 

0 


denotes time average value; also, particular index 

denotes oscillatory quantity 

denotes vector quantity 

denotes i^^ iteration 

denotes amplitude of complex quantity 

denotes steady-state value; also, source coordinate 

denotes chamber condition 

denotes droplet condition 

denotes droplet condition 

denotes fuel 

denotes gas 

axial mode index 

m,il,q 

circumferential mode index 
radial mode index 
denotes oxidizer 
denotes radial components 
denotes spray parameter 
denotes throat condition 
denotes axial component 
denotes angular component 
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INTRODUCTION AND SUMMARY 


As requirements for new rocket engines arise, selections of the propellant combi- 
nation and general hardware configuration to be used are normally required. The 
selection processes should include an evaluation of the likelihood of encounter- 
ing combustion instability problems during engine development and the expected 
difficulty of obtaining adequate stability through the use of instability suppres- 
sion devices, such as acoustic cavities. In addition, such factors as the appli- 
cability of particular stabilization devices, their size requirements, and special 
difficulties with particular propellants must be considered. 

The purpose of the investigation described herein was to analytically predict the 
relative stability of various propellant combinations of interest to NASA-JSC when 
used with hardware configurations representative of current design practice and 
with or without acoustic cavities to improve stability. Originally, consideration 
was given to propellant combinations of the LOX/hydrocarbon , LOX/amine, and the 
NTO/amine families. However, because of subsequent diminished interest in the 
former two families, attention was later concentrated on the NTO/amine family 
after the program was under way. The investigation included a literature search, 
to define the known stability- related characteristics, and stability analyses, 
employing two combustion instability models and a model to describe the effects 
of the acoustic cavities. The evaluation was done in a manner that was not spe- 
cific to any particular engine or hardware design but which provides general 
information necessary to rationally incltide stability factors in the propellant 
selection processes for a range of engine applications. In addition, preliminary 
development was completed of a generalized acoustic cavity model that includes 
steady flow and combustion effects not considered in the previously developed model. 

Early in the program, a thorough, but not exhaustive, review of the literature was 
made, to locate and summarize available information relating to stability. Infor- 
mation concerning the propellant combinations of interest and, also, similar pro- 
pellant combinations, along with any special characteristics of the propellants 
(such as a tendency toward "popping'' or to form an explosive adduct) , which may 
affect stability, was sought. The available stability-related data were summar- 
ized and related to injector and engine characteristics as well as possible. 

Although considerable information was obtained, the results suggest a complex 
interaction between important effects, Probably the results may be best used with 
the aid of an analytical model to permit separation of diverse effects. 

Two combustion instability models were used to predict the variation in engine 
stability with hardware configuration and propellant combination, exclusive of 
acoustic cavity effects. One of these was a Priem-type model; this type of insta- 
bility model has been used extensively at Rocketdyne and elsewhere for many years. 
The second model, which is applicable to liquid propellant engines but which has 
been developed for analysis of turbojet engine afterburners, has been called the 
Northern Research and Engineering Corporation (NREC) model (Ref. 1 ) . The NREC 
model is based on an analysis of liquid propellant combustion instability by 
Culick (Ref. 2), In addition, a model for the steady-state combustion was used 
to provide needed input for these stability models. Results from these stability 
analyses were used to develop parametric representations, as much as possible, to 
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show the relationship o£ stability to key parameters Csuch as dropsize and burning 
rate). Parametric representations were used to avoid restricting the results to a 
particular set of in j ector/chamber conditions and dimensions. Results from the 
analyses were used to establish the relative stability of the various configura- 
tions and to establish stability trends with changes in design parameters. 

The Priem-type stability analysis was done for a relatively broad range of condi- 
tions and extensive plots have been developed of the variation of the stability 
index, Ap, with important physical parameters. These plots may be used to a assess 
the stability of new hardware configurations with minor additional calculations. 
Calculations were also made to show the effects of varying the choice of analysis 
parameters. The results show the importance, when performing stability analyses, 
of the choice of parameters to be held constant. In addition, the results suggest 
that the propellant mass flux within an engine has a very strong influence on sta- 
bility. The greatest changes in stability were predicted when (1) chamber pres- 
sure was increased by increasing the propellant flow through a fixed set of hard- 
ware (worsened stability) or (2) contraction ratio was increased with a fixed 
thrust (improved stability). 

The analysis done with the NREC stability model was less extensive than that done 
with the Priem-type model. However, the results from the NREC analysis appear to 
agree in a rough qualitative way with those from the Priem-type analysis. The 
degree of agreement appears compatible with the major differences in approach 
that have been used in development of the two models. The NREC model, or similar 
approaches, appears to be a valuable method of analysis for liquid-propellant 
engines. The Priem and NREC models complement each other because each includes 
important factors not included in the other. Further work with the NREC model is 
recommended. 

Also, preliminary development of a model was undertaken to aid in the design of 
acoustic cavities, which specifically included the effects of the nozzle, combus- 
tion, and steady flow. This model was based on a combination of the concepts used 
in the existing cavity damping model (which does not specifically allow for these 
effects) and those used in the NREC model. Preliminary development of this model 
was completed, the effects due to pressure-coupled combustion response (velocity 
coupling effects remain to be added), nozzle effects, and steady flow (as a uni- 
form approximation), as well as the acoustic cavity being included. The limited 
computational results obtained thus far show the importance of interactions between 
the effects of the cavities and those due to the nozzle, combustion, and steady 
flow, under some circumstances. Because of these interactions, the need to ade- 
quately allow for them in cavity design is evident and additional work is recommended. 
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STABILITY EXPERIENCE REVIEW 


Available information relating to combustion instability, particularly concerning 
the LOX/hydro carbon, LOX/ amine, and N204/amine and similar propellant combinations 
has been reviewed to summarize and assess available information. The principal 
effort involved review of the proceedings from combustion and combustion instabil- 
ity meetings. Although a computerized library search resulted in the acquisition 
of only a few pertinent references, the review was aided by consultations with in- 
formed Rocketdyne personnel. 

Results from this review are summarized in Table 1, which shows the combustion 
stability experience obtained with various propellant combinstiors during a number 
of engine development programs. The references quoted in Table 1 are Ref. 3 
through 37. Abbreviations used in Table 1 are defined at the end of the table. 

The preponderance of stability experience with LOX/RP-1 has been accumulated at 
Rocketdyne. This propellant combination was used in the Jupiter, Thor, Atlas, 
Saturn IB, and Saturn V engine systems. The earliest engine development programs 
did not inlcude dynamic stability testing (i.e., no stability rating devices were 
used) but significant testing was done to eliminate spontaneous instabilities. 

The MB-3 engine, for the Thor vehicle, required an injector face baffle to elimi- 
nate spontaneous instabilities. A baffle was also employed in the MA-5 engine 
used for the Atlas booster (which was essentially identical to the MB-3); also, 
this was the first LOX/RP-1 engine for which dynamic stability was verified. An 
identical baffle configuration was employed in the H-1 engine to achieve dynamic 
stability. Much difficulty was encountered in the development of a dynamically 
stable F-1 engine. Among these difficulties were the elimination of low-frequency 
(nonacoustic) instability modes and a high sensitivity of acoustic modes of in- 
stability to minor changes in the injector face orifice pattern. In the final 
version of the F-1 engine, a 13- compartment, 3-inch-long baffle was used on the 
injector face. The F-1 engine is relatively unique compared to all other engines 
listed in Table 1, however, because of its large chamber diameter (39.2 inches) 
and its small contraction ratio (1.25). No more recent stability vork was found 
concerning engines in which the LOX/RP-1 propellant combination was used. 

With the exception of a very limited number of tests with LOX/UDMH in a high- 
pressure, two-dimensional model of the F-1, no experience with LOX/amines was 
found . 

The NTO/ amine propellant combinations have been employed extensively in recent 
years. The majority of applications have been in engines of moderate size (-10- 
inch diameter) . When these propellants have been employed in chambers with diam- 
eter greater than 10 inches, dynamic instability problems have invariably occurred 
These problems have been solved in a number of cases (GEMSIP, Apollo, SPS, etc.) 
through the use of injector face baffles. In other cases, such as the 7.79-inch- 
diameter Lunar Module Ascent engine, acoustic absorbers have been used in addition 
to (or in place of) injector face baffles to achieve dynamic stability. Moreover, 
no maximum engine size is evident below which no instability is encountered. For 
example, instability problems were encountered in very small combustion chambers, 
viz., the 3- inch-diameter Rocketdyne RS-14 (9000 Hz instability) and the 1.6-inch- 
diameter Thiokol C-1 (17,000 Hz instability). 
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TABLE 1. SUMMARY OF RESULTS FROM STABILITY EXPERIENCE REVIEW 
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3 

None 




Bomb 




Stable 

9 









LOX 

UDMH 

T/LD 

ID 

17 0 

18 F 



-30K 

-360 

2. 1-2.6 

1100 







2D 

20x1.5 

40 

1.25 


Yes 



None 
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Concluded RP-1 more stable than UDMH, particularly in the LDX-trlplet injector^ Rf-1 flow 
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damped UDMH instability in one case 



LP2D 



LOX 

RP-1 

LD 

LD 


0.073 

0.067 

5K 


2.3 

425 







2D 

21.5x1 

24 






None 




Press 


Rtflhio 


10 















i 




















Pulse 





( For uniplanar injector with short fuel impingement distance, NTO/UDMH spontaneously unstable. 



LP2D 



NTO 

UDMH 

LD 

LD 


0.073 

0.067 1 

3K 



' 42 s 







20 

21 . 5x1 

24 






None 




Press 




10 

(For biplanar injector with long fuel imp,, iP rating was: IjDX/RP- 1 = 0.92 and NTO/UDMH - 1.00 



































Pulse 





) Concluded NTO/UDMH no more tractable than LOX/RP-1 



Barrel 



LOX 1 

RP-1 

LD 

LD 


0.05S 

0.055 

12K 


2,3 

320 ! 

255 

155 





Cyl 

21 

6 

24 

18.3 

f 





None 
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1 






















I dP Rating: IXlX/RP-l = 1,00 Conclude: NTO/UDMH has lower inherent 








NTO 

UDMH 

LD 

LD 






2.14 

320 1 

200 

200 





Cyl 



1 






None 

1 










> NTO/UDMH = 0.45 resistance to instability 


























1 
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J NTO/UDMH = 0.13 than LOX/RP-1 








NTO 

UDMH 

LD 

LD 






2.55 

320 

255 

155 





Cyi 



T 

T 





None 

i 
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Horii. 

Teat Stand 

LOX 

UDMH- SO 






USX 

! 















None 



None 




None 


See N 

ote 


4 

\ 



Horiz. 

Test Stand 

LOX 

25 UUMH 
75 N 2 H 4 






^ USX 
















None 



None 




None 


Unstable 


4 

(shutdown triggered instability. Concluded N.H. either with UJX or NTO is far from tractable. 








LOX 




































/UDMH felt to bo considerably better. 

) 








^2^4 






U5K 
















None 



None 




None 
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NTO 




1 



USX 
















None 



None 




None 
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SE 

- 6 , Geitlni 
Re-entry 

NTO 

MMH 

UD UD 

Splashplate 

4 

0.0225 

0.016 

23.5 

5.8 

1.3 

> 150 

36 

85 

49 

97 



Cyl 

0 

70 


3.7 

93.6 

1 None 



None 


i 


None 


Stab 

le 


13 




SE 

- 

. 

7, Gemini GAMS 

NTO 

Mk&i 

UD UD 

Splashplate 

1 

0.024 

0.018 

79 

4.9 

1.2 

137 

21 

56 

35 

56 



Cyl 

1 

25 

3.74 

3,4 

87.6 

None 



None 




None 





13 




SE 

- 

7, GenlniOAMS 

NTO 

MKH 

UD UD 

Splashplate 

16 

0.024 

0.016 

94,5 

5.9 

1.2 

140 

30 

45 

43 

67 



Cyl 

1 

25 

4.0 

3.1 

65.9 

None 



None 




None 





13 




SE 

- 

a, Apollo RCS 

NTO 

MMH 

UD UD 

Splashplate 

16 

0.026 

0.021 

93 

5.8 

2.1 

140 

40 

57 

44 

52 



Cyl 

1 

25 

4.3 

3.1 

91,9 

None 



None 


■ 


None 
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SE 

-5, Titan 
Trans tage 

NTO 

j A-SO 

UD UD 

Splashplate 

4 

o;o24 

0.021 

25 

6.25 1 

1.56 

150 

31 

t 

33 

46 

60 



Cyl 

0 

70 

2.77 

4.45 


None 



None 




None 
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SE 


5, Titan fi 

NTO 

A- SO 

UD 

UD 

8 

0.024 

0.020 

45 

5,63 

1.56 

140 ^ 

2S 

33 

41 

62 



Cyl 

1 

06 

2.39 

5 

93.2 

None 



None 




None 





13 




Transtage 




1 





















1 















f 

SE 


10, 

LM Descent 

NTO 

A-50 

F-O-F 

Unlike Triplet 

165 

0.067 

D.Q4S 

! 10. 5K 

j 

63.6 

1.6 

ISO 

40 

40 

46 

57 

T^-40 to 
T^.lOO 


Cyl 

11.35 

15 

2.5 

97.2 

© 

1.75 

3 

None 


i 1 

Pops 

Bomb 


Marginal^ 

Marginal** 

14 

*Pop5 triggered instability which damped 
••"Almost all" damped 



RS 

-18, 

LM Ascent 

NTO 

A- 50 

F-O-F 

Unlike Triplet 

96 

0.0544 

0.0326 

■ 3.SK 

j-36.5 

1.6 

120 

51 

29 

46 

5? 


Cyl 

7 

79 

11 

2.9 

96.0 

© 

1.75 

3 

QWR 

0.6 L X 
O.T W 

‘ 5.2 


Bomb 


Stable 

Stabl e 

5 

Originally designed with baffles (IT) and cavities (IR, 3T) 


foldout feame 

( 


FOLDOUT FRAME 


K>LD 0 UI IfHAMM, 



(?- 935-3 



TABLE 1 CContinued) 


CompaTiy 


Rocketdyn^ 


A«toJ Bt 


Engi:\e 


RS-13, LM Ascent 


RS-U, PBPS 


RS-21, Mars 
Mariner 


Adv. Agena Mod. 
Mod. S533 


Titan I 
Booster 


Titan II, III 

Geninl, 1st Stag* 


Titan I Suatainer 


Titan II. Ill 
Gemini j 2iid Stage 


Titan 2nd Stage 
GEMSTP 


Apollo SPS 


Transtage 21- ID 
Transtage 21-11B 
Transtage 21- UB 
S Spxid E>ulse Motor 


LOX 

NTO 


MTO 

NTO 

NTO 

NTO 


RP-l 

A-SO 


Element Type 


Quid iter 


Fviel 


Number 

of 

Primary 

Elements 


F-o-r 

Unlike Triplet 


F-O-F 

Unlike Triplet 


LD 

LU 


U) 

LD 


LD I LD 
Quad let 
Quadlet 


2-on-2 

Unlike 

Impingement 


UD 


UD 


Quadlet 

Alt Rings Quad. 
5 F-O-F Triplet 
Alt Rings Quad. 
^ F-O-F Triplet 


LD 


LD 


O-F-0 
Triplet 
Like Doublets 
Kith O-F-0 
Fan Trip 

LD I LD 
O-F-0 
Trip] et 
Like Doublets 
With O’F'O 
Fan Triplet 


S60 0 
610 F 


563 0 
516 r 


32S 0 
392 r 


.^36 

336 

336 
40 0 
IS F 
20 

7S 


40 0 
18 F 


20 

78 


Orifice Diametei , 
inch 


Oxidizer 


0.0SO4 0 
D.0362 I 


0.0270 0 
0.0250 I 


0,119 

O.llD 

0.085 


0.036 

0.036 Q 
0.051 T 
0.036 Q 
0.051 T 

0.078S 

0.104 

0,028 


0.07B5 
0. 104 
0.028 


0.0397 0 
0.0319 1 


0.0175 0 
0.0205 I 


0.037 


0.041-0.077 

0.0292 

0.0292 

0,0292 

0.0785 

0.104 

0.035 

D.0785 

0.104 

0.035 


Thrust, 

pounds 


SOK 

lOOK 


8.15K 

8.1SK 

8.15K 

4.5K 


Thinjst 
pe r 

Element , 
Ibf 


111 

100 


200 

200 

30 

40 

40 

40 

24 

24 

24 


Mixture 

Ratio 


1.57 

1.6 


2.25 

2.0 


Chamber 

Pressure, 

psia 


1 . 6 * 

2.0 

2.0 

2.0 

0.9-1.44 


120 


682 

827 


Injection 
Pressure 
Drop, psi 


100 

100 

100 

300 


Inj ection 
Velocity, , 
ft /sec 


-40 

w40 

-40 


fou>out frame 

I 


;'^4o 

~45 

.-45 


Injection 

Temperature, 

F 


Ta- 40 to 
Ta*^100 


31-73 

24-85 

24-85 


32-60 

57-86 

57-86 


General 

Chamber 

Configuration 


Cyl 


Cyi 


Cyl 


Cy] 


Chamber 

Diameter, 

inch 


21.6 

21.8 

14.2 

14.5 


17.6 


11.91 

11.91 
11.91 
12.8 



Injector 
to Throat 
Distance , 
inch 


17.72 

17.72 
17,72 


Contraction 
Ratio 


2.53 

2.53 

2.53 
10.3 


97.8 
97,2 

98.9 


© 


None 
Tabs 
0.5 K 
Tabs 
0.5 K 

None 


See 

Note* 


None 

© 

None 

© 

None 

None 

None 

None 


None 

© 

None 


@ 

© 

© 

None 

None 


Baffle 

Length, 

inch 


1,75 
1 .75 
1.0 
0.75 

0.6 L 
0.6 L 


0.4 to 0.4 

1.25 Y 
0.5 0 


1-1/4-2 


4 

4 

4* 

3.2 

3.2 


Number of 
Blades or 
Compartments 


See Note* 
5 


8 « Hub 
6 + Hub 
7 


6 Unsym. 
+ Hub 
5 * Hub 

5 + Hub 
4 


Acoustic 

Absorber 


I Absorber 
Size 


QWR 


None 

None 


Groove 
in Wall 


600 Orlf 
Liner 


None 

None 

None 

None 

None 

None 

None 


None 

None 

Yes 

None 


Absorber 
Fractional 
Dpen Area, 
pcrc ent 


0.6 L X 
0-1 W 


O.SS L X 
0.06 W 


1/16 X 

7/32 


5.2 


Pops, 

etc. 


Ign. 

Spikes 


Ign. 

Spikes 


Ign. 

Spikes 


Pops* 


Pops* 


5. 3/6. 5 


Pops* 

Pops 

Pops 






Undisturbed 

Dynamic 



bil i ty 


Stability 

Stability 



sting 

ivice 

Instability 

Mode 

(No Rating 
Dotficej 

(With Rating 
Device] 

References 

Stabi^ty Comments and Conclusions 

omb 


Stable 

Stabl e 

15 

) Originally designed with baffles (IT] and cavities (IR, 3T) 


IR, 3T 

Stable 

Unstable 

IS 

) 




Stable 

15 

\ Determined marginal baffle length 




Unstable 

IS 

) 


IT 

Stable 

Stable 

15, 16 

Could Stabilise with cavity alone; 16 percent open area marginal 

►ulse 

lode 

IT 


Unstable 

17 

1 to 2 percent Incidence rate of instable (3 percent without tabs] 




Stable 

17 

Completely stable 




Stable 

17 

Completely stable 

lone 



Stable* 

5 

•Assumed stable 


IT 

Stable 

Unstable 

18, 19 





Usually 

Unstable 

IS, 19 

*CD©®^ Not all blades extended to wall; especially unstable with wall 

^ gap. “Pops associated with propellant accumulation under baffle 


IR, 3T 


Unstable 

18, 19 

Concluded ineffective damping of IR which triggers 3T 




? 

18, 19 

Investigated but dropped 




Stable 

18,19,20 

Final injector configuration included groove in wall at chamber periphery 
to eliminate 5T 


IT 


Unstable 

Marginal 

21 

21 

1.25 inch long baffle dyn. unstable; l.S inch and 1.75 inch baffles resurge; 2 inch 


2T 


long baffle dyn, stable but long damp time. ‘Associated with accumulation under baffle 

Jomb 



Unstable 

21 

Discontinued effort 





21, 22 

Recommended but program terminated before testing 



Stable 

Unstable 

2.3 

24, 25 

Statistically stable spwntaneously 

Statistically stable spontaneously# 0.6 percent incidence of instability (during 

Bomb 

IT 

Stable 

Unstable 

operation in low MR, low Pj. region) 



Stable 

Unstable 

24 

Statistically stable spontaneously 

?.G. 

IT, 2T 

Marginal* 

Unstable 

25 

*1.3 percent incidence of instability 

lomb 

IL 

Unstable* 

Unstable 

2S 

*30 percent incidence of instability 

\G. 

8amb 

3T 

Marginal 

Unstable 

25 

*1.3 percent incidence of instability 

P.B. 

lomb^ 



Stable 

25, 26 

•Marginal length 2,5 inches. Ifumped injection at tiid-radius 




Unstable 

26 

/ 


i 


Stable 

26 

•Marginal length 2.5 inches. Flat distribution; selected concept 

Bomb 

2T 

Stable** 

Unstable 

25 

•Less than 1 percent incidence of instability. **Stati5tically stable. Injector 


Stable 

Stable 

2S 

has varying impingement distances across face 
. Mixed impingement 




Stable 

Stable 

25 

Mixed impingement; selected concept 




Stable* 

27 

•Shift in MR from 2.0 to 1.6 resulted in marked increvrierits in incidence of random popping; 
dyn. stab, was restored by counterboring ox, orifices and increas.ing baffle length 


2T 

Unstable 


25, 27 

Pops triggered sustained instability 



Stable 

Marginal 

25, 27 

Flight version 




Unstable/ 

Stable 

5, 24 

6.5 percent open area of absorber required for dyn, stability 

P-G. 


Unstable 

Unstable 

28 




Unstable 


28 

j 



Unstable 


28 

1 ‘Unstable at low MR. 





28 

V Concluded MMH more stable than N2H^. 




Stable 

1 Instability increases with decreasing MR for N2H4. 

1 Experimental results agreed with predicted stability trends 




Unstable 

Unstable 

26 

1 using Priem inode 1 



Unstable 


28 
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TABLE 1 (Concluded) 


Company 


Thiokol 

UTC 

Marquardt 


Engine 


8 Spud Pulse 
Motor 


Transtagc 21-lU 


Transtage 21- ill 


KC-1 

C-1 


S-IV-B Ullage 


"Zot'' Study 


UDMH-50 

^WIi 

MMH-50 

VDm~50 

MMil 

MMH-SO 

MMH 

UDMH 

UDHH-50 

MMH 

UDMH- 50 


UDMH- SO 
MMH 
UlAfH 


Element Type 


Quad let 


Alt, Rings Quadlet 
and F-O-F Triplet 


Quadlet 

UD I UD 


WUmber 

of 

Primary 

Elements 


Orifice Diameter, 
inch 


0.042 
to 0.12 


0.056 Q 
0.051 T 


0.039 
to 0.113 


Thrust, 

pounds 


IDO 

8K 

1 .75K 


Thrust 

per 

Element , 
Ibf 


13S 

10 

15 

IB 


Mixture 

Ratio 


2.0 

1.6 

2.0 

1,62 


ChaDiber 

Pressure, 

psia 


300 

400 

600 


15 to 
450 


100 

96 

100 

100 


Injection 
Pressure 
Drop, .p si 


Oxidizer J Fuel 


Inj ection 
Velocity, 
ft/sec 


40 

to 130 


30 

to 120 


Injection 

Temperature, 

F 


104-126 104-127 


General 

Chamber 

Configuration 


Cyl 


Cyl 


Cone 

Cyl 


Chamber 

Diameter, 

inch 


18 

1.60 

10 

4.75 


Injector 
to Throat 
Distance, 
inch 


Contraction 
Ratio 


U (Pc5 


3.38 

1.81 

1.84 

1.84 

2.75 


-9,5.0 

-97.8 

-91.8 

-98 

-98 


None 

© 


© 


None 

None 

None 

Added 

None 


Baffle 

Length, 

inch 


Number of 
Blades or 
Compartments 


5 Rad Tabs 
10 Rad Tabs 


Acoustic 

Absorber 


Helmholtz 
Reson . 


Absorber 

Size 


Absormr 
Hractiepi 
Open 
percei 


i Ai{c: 


lit 


Pops, 

etc. 


Pops 


Less 

Pops 


Pops 

Pops 

No 

Pops 


Stab i 1 i cy 
Rating 
Device 


Bomb 

P.G. 


Bomb 

P.G. 


Pulse 

Made 


Instability 

Mode 


IT 

IT 

4T, 5T 


Undisturbed 
Stability 
(No Rating 
Device) 


Stable 


Stable 

Stable 

Unstable 


Dynamic 
Stability 
(With Rating 
Device) 


Unstable 


Unstable 

Unstable 


Stable 

Stable 

Unstable 


Unstable/ 

Stable 


29 

32 

33 

34 

35 


36, 37 


Stability Comraents and Conclusions 


Stability increases with: Inc. P^, dec. W'j-, Inc. 

trend agrees with Priem theory 

pr^ictions 

S to nothing; unstable to pops 1/2 

S to 10 gr. P.G. 2S percent of time (1/4); unstabli 
to pops 2/7 

S to 10 gr. P.B., 15 gr. P.G, 1/7; unstable to 
pops 7/19 

Always unstable to pops 6/6 

Concluded: Inc. S with inc. MR. MMH stability 

— UDMH-50 stability until MR > 2 

— 20 gr. P.G, threshold 


20 gr, P.G, threshold 
Unstable to pops 2/3 


No correlation be- 
tween pops and pro- 
pellant temp., per- 
cent P.C, baffle cracks 


Pop study; Observed that pops inc. with inc. 

DoTif, dec. Vinj> dec. P^, dec. Tinj 

Popwise N2H4 —MMH with UDMH much better (almost no 
pops) 

Req. baffle length between 2.4 and 2.9 for dyn, 
stability 

10 percent instability occurrence rate without absorber 

Unique in that 1, 2, 5T modes were triggered but 
damped 


Unstable with S tabs, stable with 10 tabs 
Achieved stability by inc. CR 


Loudest Explosion. 2ots occur when cell press 
>1/2 VP of N2H4 at head end temperature. 

Zots aggravated by chamber wall temperature 
>head end temperature 


ABBREVIATIONS 


Ann 


annular chamber 

P.G. 

CR 


contraction ratio 

QNR 

enpt 

- 

baffle compartment 

RP-1 

Cyl 


cylindrical chamber 

S 

D 

= 

mass median droplet size 

Spent 

dyn. 


dynamically 

^A 

elim. 

= 

eliminate 

T/LD 

F 

- 

fuel 

UD 

GEMS IP 


Gemini Stability Improvement Program (Aerojet) 

UDMH 

Horit. 

= 

horizontal 

UDMH- 50 

jpl 


Jet Propulsion Laboratory 

Unsym. 

L 

■ 

acoustic cavity depth 

US 

LD 


like doublet injection element 

UTC 

LH 


lunar module 

Var 

LOX 


liquid oxygen 


L. Trip 

- 

like propellant triplet element 

W 

MMH 


monomet hyl hydrazine 


NTO 

- 

«2°C 

ll 

0 


oxidizer 

IT 

Orif, 


orifice 

2D 

PBPS 

= 

post boost propulsion system 

2T 

P 

c 

“ 

chamber pressure 

3T 


AP 
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MJX)UT FEAlffi: 


= pulse gun stability rating device 
= quarterwavc acoustic resonator, acoustic cavity 
= kerosene-type fuel 

- stable 

= spontaneous 
■ ambient temperature 

= mixed triplet and like-doublet injection elements 
= unlike doublet injection elements 
= unsymmetrical diraethylhydrazine 
= SO percent and 50 percent UDMH 
■■ unsymmetrical 
= unstable 

= United Technology Center 
' various 

» injection velocity 
= width of acoustic cavity 

- total propellant flowrate I 

- first radial mode of instability 

« first tangential mode of instability 
= two-dimensional combustion chamber 

- second tangential mode of instability 
= third tangential mode of instability 
= pressure pulse amplitude 

I 

I 




f 
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As might be expected, little information was found that would allow a direct com- 
parison between LOX/RP-1 systems and NTO/ amine systems. However, a limited num- 
ber of directly comparable tests was made with LOX/RP-1 and NTO/UDMH propellant 
combinations in Rocketdyne engines. Both propellant combinations were operated 
successfully in the Atlas sustainer, Thor/Jupiter, and E-1 engine systems but 
without attempting to determine the relative dynamic stability. In addition, both 
propellant combinations were used in a low-pressure, two-dimensional model of the 
F-1 and a so-called "barrel" chamber with inert gas pulse rating techniques being 
employed. The barrel chamber was a model of a relatively high thrust engine in 
which the chamber diameter was modelled but a high contraction ratio was used to 
reduce the thrust to a convenient level. The results indicated essentially no 
difference of stability in the two-dimensional engine but an indication of greater 
stability with the LOX/RP-1 combination with the barrel chamber was found. In 
addition, the LOX/UDMH propellant combination was found to be less stable than 
the LOX/RP-1 combination in tests with a high-pressure, two-dimensional model of 
the F-1 chamber. It should be noted, however, that in most cases no special effort 
was made to operate with the alternate propellant combination at its optimum mix- 
ture ratio and in many cases the exact mixture ratio achieved is not readily 
available. 

The NTO/SO-SO propellant combination was substituted for LOX/RP-1 when the Titan-I 
was uprated slightly to become the Titan-I I. Both the Titan-I and Titan-II were 
only "statistically" stable (i.e., vrithout disturbing the engine with a rating 
device, large numbers of tests could be made without the occurrence of spontaneous 
instabilities) . These engines were unstable to artificial disturbances introduced 
with dynamic rating methods. No inherent stability advantage of either propellant 
combination appeared evident. The desire to achieve a dynamically stable Titan 
second-stage engine led to the GEMSIP (Gemini Stability Improvement Program) Pro- 
gram which culminated in the achievement of dynamic stability through an increase 
in thrust per element and the addition of injector face baffles. 

The AFRPL performed a comparative stability study of the Transtage engine with 
N2H41 N2H4/UDMH (50-50), MMH, and N2H4/MMH (50-50) fuels and NTO oxidizer. With 
pulse guns being used for stability rating, the stability with the N2H4/UDMH (50- 
50) and MMH fuels was approximately equal and was significantly better than that 
obtained with either the N2H4 or N2H4/MMH (50-50) fuels (which were also conducive 
to "pops") . 

The susceptibility of amine-type fuels to pops, ignition spikes, etc., has often 
been of concern. While pops have been observed with LOX/RP-1 propellants, these 
pops have been attributed to entrainment of air in the propellant stream (Ref. 5). 
Pops observed with NTO/amine propellants have generally been attributed to blow- 
apart or stream separation (Ref. 51). Pops with NTO/amine propellants often re- 
sult in sustained instability and/or hardware damage. A recent study of popping 
is reported in Ref. 30 in which the popping tendencies of N2H4, MMH, and UDMH were 
compared under varying injection conditions; UDMH was found to have less popping 
tendency than N2H4 or MMH, which were comparable to each other. 

Minton and Swick (Ref. 36) report the investigation of manifold explosions (which 
they called "Zots") with NTO/amine propellants. Their experiments indicated that 
the condensation of residual fuel vapor in a cold, empty oxidizer manifold followed 
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by the contracting of this detonable fuel condensate with fresh oxidizer leads to 
a manifold explosion or "Omzof'. In Ref. , other types of Zots (Chizot, Pizot) 
are described, some of which would be called pops by other investigators. Zots 
were concluded to occur only under conditions where local pressure is greater than 
one-half the vapor pressure of N2H4 at the local temperature. The .occurrence of 
Zots was found to be aggravated by high chamber wall temperatures. . Artificial 
explosions were created by adding NTO to identical amounts of various fuels. The 
"loudness” of the resulting detonations was measured and the resultant ranking may 
be an indication of susceptibility to occurrence of damaging Zots as well (Table 1] 

Xhe literature review indicated that UDMH is the least likely of the amine family 
to exhibit detrimental popping. MMH and (N2H4/UDMH, 50-50) are about equally sus- 
ceptible to detrimental popping and are ranked as worse than UDMH but better than 
N2H4/MMH (50-50) and straight N2H4 which are considered quite susceptible. 

Of particular interest to this program is the work of Abbe et al. at AFRPL (Table 
1)) during which the Priem stability model was used successfully to: (1) predict 

the relative stability of an NTO/MMH engine design with an NTO/N2H4-UDMH (50-50) 
engine design (Ref. 27) and (2) to investigate the effect on stability of various 
design and operating conditions (Ref. 28). Abbe concluded that interactions be- 
tween the various processes occurring in a combustor are so complex that an ana- 
lytical model is necessary to predict the overall effect of a change in even a 
single design or operating parameter. 


R-9353 

8 



PRIEM MODEL INSTABILITY ANALYSIS 

DESCRIPTION OF THE MODEL 


The Priera-type instability model has been used extensively to predict engine sta- 
bility. This model is based on numerical solution of a system of equations describ- 
ing the combustion/ flow field within a small annular region chosen as representative 
of the combustion chamber. The model considers basically only one dimension and, 
therefore, can approximate tangential modes but not radial, longitudinal, or coupled 
modes. Although this model has some limitations, it does include the detailed spray 
combustion processes in a quasi-steady sense and it does account for nonlinearities. 
Numerical solution of the equations gives a nondimens ional overpressure, - Ap/pc, 
required to initiate a high-frequency instability (transverse mode) as a function 
of several nondiraensional parameters. This critical overpressure, Ap, is deter- 
mined from a series of computer experiments in which the oscillatoiy combustion 
response (history), subsequent to various initial disturbance levels, is calculated. 
The disturbance that will cause an oscillation that neither grows nor decays is the 
critical overpressure, the Ap. Conceptually, this is similar to bombing an engine 
with varying bomb sizes until the minimum size is found that will cause a sustained 
oscillation. 


According to the Priem model, the stability index Ap is a function of several 
parameters : 


Burning Rate Parameter, ^ 

Velocity Difference, AV* = I'^gas “ drops 

Mass Accumulation Parameter, MAP = M /w. . t 

vap inj w 


10 


-1/3 /D^^cp\-l/2 


Small AV* Parameter, K, = ■=— Sc 

13 0 

' / o 

/3 Cp R M \ 

Droplet Drag Parameters, K 2 = ^ ' 4 ' / ' pDj j , and 


K, = 


'3 Cp Rp 


3 \4 (pD)^^ 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


where the subscript o refers to steady-state (stable) conditions. The values of 
these parameters may be calculated from the results of steady-state combustion 
model calculations. Generally (for most propellants and injector configurations), 
the most influential parameters are the burning rate parameter, and AV', the 
axial velocity difference. 


A generalized plot of the calculated relationship between Ap, ^ and AV' is shown 
in Fig. 1. The parameter AV varies with axial position in the engine and usually 
has a minimum value near the injector. The model predicts that this region, cor- 
responding to a minimum AV, is the most unstable or most sensitive zone and, 
therefore, the analysis is generally done for this location only. Further, although 
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the value of AV may be predicted to be zero at this location, it is generally 
agreed that a minimum value of 0,01 is more reasonable because of turbulence ef- 
fects and nonuniformity of dropsizes. Thus, the stability analysis is commonly 
made with AV* = 0.01 and for the sensitive zone. 

The relationship between Ap, and Ki, for AV* = 0.01 is shown in Fig. 2. The 
parameter Ki arises through consideration of a more exact burning rate expression, 
than was originally used by Priem, for small values of AV (Ref. 38). 

The relationship between Ap, SB, and MAP is shown in Fig. 3. The effect of the 
MAP parameter is to raise the stability limit at locations where a large propor- 
tion of the propellant has already been burned, and a small amount of unvaporized 
propellant is left to sustain an instability by burning. The most sensitive zone 
of instability (minimum AV*) normally occurs at an axial location where a large 
proportion of the propellant is unburned (large MAP); therefore, the MAP parameter 
is generally of little importance. 

The relationship between Ap, SB, K2 and K3 is shown in Fig. 4. These parameters 
arise from a modification to the Priem model to include the effect of droplet drag 
(Ref. 38). These curves show that Ki and K2 can significantly change Ap, indicat- 
ing the in^ortance of droplet drag. 

Results from a steady-state combustion model are required for the Priem-type anal- 
ysis. Such steady-state models have been used extensively and several model varia- 
tions are available. Generally, in these models, the mass and mixture ratio dis- 
tributions created by the injector are assumed uniform to avoid the complexities 
of handling the nonuniform case. Moreover, in most cases it is satisfactory to 
employ a model based on an evaporation coefficient rather than including droplet 
heating. This approximation is good as long as the chamber pressure is signifi- 
cantly below the critical pressure of the propellants. 

Steady-state combustion models are based on numerical solution of the differential 
equations describing the spray combustion processes. These equations apply to 
arrays of discrete droplets and, generally, are restricted to one-dimensional flow. 
Therefore, the calculation is usually begun at a position downstream from the in- 
jector where the equations are expected to apply. The calculation then proceeds 
downstream in a stepwise fashion to the nozzle throat. Iteration procedures are 
employed to ensure that the flow satisfies sonic conditions at the throat. The 
input required by the steady-state models includes: 

• Chamber geometry (diameter, length, contraction ratio) 

• Combustion gas properties (p, T, p, k, etc.) as a function of mixture 
ratio (obtained from equilibrium performance calculations) 

• Axial position where calculation is to begin 

• Initial condition of spray and gas at starting position 

• Spray dropsize and dropsize distribution 

• Spray velocity 

• Spray temperature 
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Figure 2. Priem-Type Stability Map Showing Effect of Parameter 
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Figure 3. Effect of MAP Parameter on Predicted Stability Limits 
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BURNING RATE PARAMETER, SB 


Figure 4. Effect of K2 and Kj Parameters on Predicted Stability Limits 
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• Gas composition (MR) 

• Gas flowrate (percent of spray burned) 

• Gas pressure 

• Initial dropsize 

• Initial drop velocity 

• Gas velocity 

• Gas composition (percent of fuel burned, percent 
of oxidizer burned) 

• Initial overall percent burned (vaporized) 

Output from the steady-state combustion model includes the following variables as 
a function of axial location: 

• Dropsize 

• Drop velocity 

• Gas velocity 

• Gas composition (fraction of fuel burned, fraction of oxidizer burned) 

• Gas physical properties 

• Overall fraction burned (vaporized) 


The Priem model input parameters can be calculated from the steady-state combustion 
model output through use of the following relationships: 


m = 


Ac 

Az 


where e is the overall fraction burned. 


M 


vap^ 


= { 


il?r * 


MR+1 } ^o 


(7) 


( 8 ) 


where Xp are the burned fractions of oxidizer and fuel, respectively. Also: 


w. . 

■ =-^1" 

t = ^15. 

W ~ C 


(9) 

( 10 ) 




27.0 Re 
0.271 Re' 

2.0 


•0.84 

0.217 


Re < 80 
80 < Re < 10^ 
Re >10^ 


( 11 ) 
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where Re = 


Adrl^ 


P 


The use of the Priem-type model in conjunction with a steady-state combustion 
model allows a relatively thorough evaluation of variations in engine stability 
with the parameters of interest. These parameters can include dependent param- 
eters such as dropsize as well as operating and design variables such as chamber 
pressure or contraction ratio. 

STEADY-STATE COMBUSTION MODEL CALCULATIONS 


In formulating the analysis scheme, particular care was taken to group the effects 
of the many independent variables such that a parametric analysis would not require 
excessive computer time or manual labor. Moreover, the results must be expressible 
in a reasonable number of graphs, tables, or example problems. A summary of this 
analysis scheme and selection of analysis parameters is described below. 


A so-called k’ or evaporation coefficient steady-state model was used with the 
following input parameters: 

• Chamber diameter 

• Chamber contraction ratio (throat diameter), CR 

• Chamber length 

• Chamber pressure, p^ 

• Mixture ratio, MR 

• Initial spray velocity, 

• Initial spray dropsize, D- and D /D^ 

I OX X 

While the chamber and throat diameters are actually input into the steady-state 
model, the model is one-dimensional and actually considers only the total propel- 
lant mass flux (flowrate per chamber cross-section area) w^/A^. Also: 


“t PA4C Pc*c 

A ■ C*A ■ c*CR 
c c 


( 12 ) 


This expression was used to eliminate the chamber diameter from the list of key 
parameters shown above. 

Further, propellant dropsize was treated as an independent variable, which effec- 
tively separates some of the effects of injector orifice, injection pressure drop, 
liquid propellant properties, and chamber contraction ratio. 

The coding of the k’ computer model was modified to directly calculate the follow- 
ing Priem parameters from the steady- state model results: m, Myap^* ^injo’ 
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and K^. In addition, the following Priem- related parameters were calculated by 
the model as functions of axial position: 

^/R 

MAP(R) 

K2/R 

K 3 /R 

The effect of chamber size (radius, R) on the Priem parameters and thus, the sta- 
bility limit (Ap) , was calculated through use of Fig. 2 without additional 
steady-state model results being required. 

Additional practices were followed to reduce the number of variable steady-state 
model input parameters from the nine listed earlier t£ a ^otal of six, viz., pro- 
pellant combination, CR, pc, Vinj (or Apinj), Df and Dox/Df. The chamber length 
was held constant at 9 inches (injector face to beginning of nozzle convergence). 
The combustion efficiency (he*) was assumed to be 96 percent. Mixture ratio was 
fixed for each propellant combination (thus fixing the theoretical c*); a mixture 
ratio of 1.60 was used for both NTO/MMH and NTO/50-50 propellant combinations. 

The assumed initial percent burned at the computational starting position was 1.0 
percent of each propellant. The spray velocity was calculated from selected values 
of with Cd = 0.89 being used. A specialized form of the Nukiyama-Tanasawa 

dropsize distribution was used which allowed the distribution to be calculated 
from a specific mass median dropsize. This distribution is based on work by Ingebo 
as described in Ref. 40. The distribution function is: 

(4.63/D)^ (D^/120) exp (-4.63 D/D) D < 2.00 D 
0.0 D > 2.00 D 

where D = or Df and V is the volume fraction of the spray. An axial step 
size of 0.020 inch was used in the computer calculation procedure. The inclu- 
sion of droplet heating is the k' combustion model was not attempted. 

A total of 240 computer runs with the steady-state combustion model was made for 
the NTO/MMH and NTO/50-50 propellant combinations. These computer runs encom- 
passed the entire desired matrix of steady-state combustion conditions to be in- 
vestigated for those propellant combinations. The steady-state model runs for 
each propellant combination included all possible combinations of the following 
variables: 

Chamber Pressure, p^. = 125 psia, 200 psia 

Injector Pressure Drop, Ap = 0.2 p^., 0.4 

Combustor Contraction Ratio, CR =2. 0:1, 3.0:1 

Fuel Dropsize (Mass Median), Df = 60, 100, 200, 300, 400 microns 

Ratio of Oxidizer to Fuel Dropsize, D^j^/Ff = 0.5, 0.75, 1 .0 


dD 
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RESULTS FROM PRIEM-TYPE ANALYSIS 

Steady-state combustion parameters needed as input to the Priem instability model 
were tabulated from the steady-state model results as functions of distance and 
chamber radius, if appropriate, for each case. The location_of ;^he most sensitive 
zone was found to shift further downstream as the variables Df, Dq^» CR, Apinj and 
pg were increased. For operation at D£ = 60 microns, Dqx ~ 30 microns, CR = 2, 

Api ^ = 0.20 Pc and pc = 125 psia, the most sensitive zone was located about 0.02 
incn downstream of the starting point for the computer calculation^. The maximum 
shift in the location of the sensitive zone was found to occur at Df = 400 microns, 
Dqx = microns, CR = 3, Apinj = 0*40 pc and p^ = 200 psia. The sensitive zone 
was located at slightly more than 2.00 inches from the starting point for these 
conditions. Although the value of AV* at the most sensitive zone was assumed to 
be its "turbulence," or minimum, value of 0.01, the actual numerical value was 
usually only 0.002 to 0.004. 

In all cases, the numerical values of MAP (representing the effect of remaining 
unvaporized propellant) and K2 and K3 (representing the effects of droplet drag) 
were found to be of such a magnitude that their influence on predicted stability 
was negligible for the entire matrix of computer runs. This is not unusual for 
storable propellants. Calculation of the Priem stability index, Ap, was thus 
simplified considerably, being a function of only the non-dimensional burning- 
rate parameter, the parameter Kj (the correction for small AV*) and the non- 
dimensional AV , which was assumed to have a turbulence level value of 0.01 at 
the sensitive zone. 

Two methods of calculating the burning rate parameter were used. One method in- 
volved the use of the local value of the fraction burned per unit length (m) at 
the most sensitive zone. With the second method, the average value of m over the 
first 2 inches of the combustion zone was calculated, which included the "most 
sensitive" zone for all computer cases and which would avoid high local "point" 
values whose magnitude is sensitive to poorly defined steady-state model input 
parameters such as initial percent of propell'ant burned. The local value of m 
at the sensitive zone was almost a factor of 10 larger than the averaged value of 
m over the first 2 inches of chamber length for small dropsizes while for large 
dropsizes it was roughly one-half the averaged value. 

The calculated value of Kj increased with decreasing p^., Df and Dqx hut was inde- 
pendent of CR and Api^j . A maximum value of = 0.204 was obtained for p^ = 125 
psia, Df = 60 microns and =_30 macrons. A minimum value of Ki = 0.041 was 
obtained for pc = 200 psia and Df - Dqx “ 400 microns . 

Values of the Priem stability index, Ap, were determined for the previously men- 
tioned ranges of chamber pressure, injection pressure drop, contraction ratio, 
and fuel and oxidizer dropsize for both the NTO/MMH and NTO/50-50 propellant com- 
binations at combustor radii of 2, 4, and 6 inches. Figure 2 was used in con- 
junction with the numerical values of and Kj obtained from the computer print- 
out to obtain the Ap value. 

The results from this analysis are summarized in Appendix B, examples of which 
are shown in Fig. 5 and 6. These results have been placed in an appendix 
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because of the large number of figures involved. Results obtained with the local 
value of the burning rate are presented in Fig. B-1 through B-16, while results 
obtained with the averaged value of the burning rate over the first 2 inches of 
axial distance are presented in Fig. B-17 through B-32. The dotted lines in the 
figures r^resent an extrapolation of the stability curve shown in Fig. 2 to a 
value of ^between 1.0 and 10. For purposes of extrapolation it was assumed that 
the stability limit remained invariant beyond ^ - 1.0, i.e., Ap > 1.0) = 
ApCi^=1.0). Attempts during previous Rocketdyne programs to extend the range of 
the stability limit plot in the direction of greater ^ were unsuccessful because 
numerical instability was encountered in the computer solution of the programmed 
equations. While Priem's original work (Ref. 41) showed a bowl-shaped curve, with 
Ap increasing with increasing ^ for ^ > 1 , more recent work at Rocketdyne indi- 
cates that it may not increase, the apparent increase being due to numerical in- 
stabilities. Therefore, the dotted lines are believed to be upper limits for the 
stability index. . 

Many of the Ap curves shown in Fig. B-1 through B-32 exhibit a minimum in Ap. 

This results from the influence of Ki, which is dominant over the influence of the 
^ . It occurs most often when the value of ^ is relatively large and the varia- 
tion of Ap with ^ is small (see Fig. 2 ). Although this minimum in Ap implies 
a "worst” dropsize, the curves on which the results are based have been extrapo- 
lated somewhat and this conclusion may be invalid. 

The Priem-type stability analysis showed no significant difference in stability 
between the NTO/MMH and NTO/50-50 propellant combinations. The effect of varying 
each of the other parameters on stability, the remaining parameters being con- 
stant, can easily be determined by comparison of Fig. B-1 through B-32. These 
results indicate that stability increases (improves) with increasing dropsize (for 
fixed CR, Apinj/Pc> Pc» Dox/Df» and chamber diameter), contraction ratio (for 
fixed dropsize, Api^j/pQ, p^. and chamber diameter) and fractional injection pres- 
sure drop as well as with decreasing chamber pressure and chamber diameter. How- 
ever, care must be exercised in making and using such comparisons because the re- 
sults depend on the parameters being held constant. The parameters being held 
constant for the results just noted imply hardware changes to allow them to be 
held constant. Different results would be expected if, e.g., chamber pressure 
were changed with fixed hardware. 

The parameters chosen for the stability calculations (chamber pressure, fractional 
injection pressure drop, contraction ratio, chamber diameter and dropsize) were 
chosen because the selection minimized the number of parameters considered and 
emphasized the important physical processes. Nevertheless, the results (Fig. B-1 
through B-32) may be used by interpolation to estimate the Priem stability index 
and, thus, relate stability to various hardware configurations and operating con- 
ditions. Such calculations have been made that illustrate the use of the curves. 
The principal additional information that is needed for this assessment is some 
method of calculating dropsizes of the propellants. 

The influence of varying certain design parameters on stability with alternate con- 
straints, relative to a reference case, was calculated. The reference case was: 

pc = 125 psia 

D =8.0 inches (chamber diameter) 
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CR 

= 2.0 

^Pinj/Pc 

= 0.40 

Df 

= 300 microns 

Dox/Df 

II 

O 

0 

cn 

MR 

= 1.6 


The local stability limit for this base case = 0.0164) is given by 

Fig, B-5. 


The parameter constraints that were used to relate the various model parameters to 
the design or operating parameters are summarized below. Thrust is related to 
chamber pressure and contraction ratio through the thrust coefficient: 


w C c* p A C 




CR 


(13) 


Injection velocity can be expressed as either a function of Apj^j^j or of flowrate 
per orifice and orifice size where the flowrate per orifice is, xn turn, function 
of thrust per element: 


V. . 
in; 



w . _ w . - 
orif _ orif 


(14) 

(15) 


where 

orif Cp n * 
I* c 


and C(Fg) is a function of the MR and element type that relates flow per element 
to flow per orifice. 

The dropsize may be related to such variables as orifice diameter, injection velo- 
city, and chamber pressure through any of a nimiber of correlations. The Dickerson 
like-doublet dropsize expression (see Appendix A) was used to determine the ef- 
fect of various parameters on dropsize. This equation gives the relationship 
D ~ - 0.852 d„jf0.568. 

Results from the calculations for alternate parameters are summarized in Table 2 
and are discussed below. If chamber pressure is varied at a constant thrust level, 
then either chamber area or chamber contraction ratio must also vary (see Eq. 13) 
but total propellant flow is fixed. Assuming the number of injector elements is 
held constant, either (1) the injector pressure drop, Ap — , may be fixed or (2) the 
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TABLE 


2 


PREDICTED EFFECT OF PARA^iETER VARIATIONS ON STABILITY 


Parameter Variation 

Pc’ 

psi 

Chamber 

Diameter, 

inch 

CR 

Ap. . , 
psi 

Ap. 7p 

^inj 

Inj ection 
Velocity, 

Thrust, 

Thrust 

Per 

Element, 

Injection 

Orifice 

Diameter, 

D/D (2) 

ref 

D (1^ 
f * 
microns 

Stability 

Index, 

A /A 
P p,ref 

Reference Case 

125 

8,0 

2.0 

50 

0.40 

1.0 

1.0 

1.0 

1.0 

300 

1.000 

Increase Chamber Pressure 












Increase Flowrate 

200 

8.0 

2.0 

128 

0.64 

1.6 

1 . 6 ^ 

1.6 

1,0 

201 

0.884 

Reduce Throat Area 

200 

8,0 

3 . 2 *^ 

ISO 

0.25 

1.0 

1.0 

1.0 

1.0 

300 

1.043 





180 

0.40 

1.26 

1.0 

1.0 

0.88 

230 

1.012 

Reduce Chamber Diameter 

200 

6.3 

2.0 

(SO 

0.25 

1.0 

1.0 

1,0 

1.0 

300 

1.049 





180 

0.40 

1.26 

1.0 

1.0 

0.88 

230 

0.988 

Increase Contraction Ratio 

i 











Increase Chamber Diameter 

125 

9.8 

3.0 

50 

0,40 

1,0 

1.0 

1.0 

1,0 

300 

1.299 

Reduce Throat Area 

187 . 5 ^^^ 

8.0 

3.0 

(SO 

0,27 

1.0 

1,0 

l.O 

1.0 

300 

1,067 





|75 

0.40 

1.0 

1.0 

1.0 

0.90 

248 

1.055 

Reduce Injection Ap 

12 S 

8,0 

2.0 

25 

0,2 

0,71 

1.0 

( 0,71 

1.0 

403 

1.140 









1 1.0 

1,19 

419 

1.189 


(1^ cases 

(2} Referred to reference case; A - - 0,0164 (local value) 

p, ret 

(3) = 0,0164 (local value); relative values greater than indicate improved, stability 

(4) Only case in which thrust was allowed to vary 

(5) Steady -state combustion model was not run for these cases; obtained by interpolation 




fractional injection pressure drop, Api^j/pc, may be fixed, with corresponding 
changes in injection velocity and orifice diameter. Results from a series of cal- 
culations for several parameter variations are simimarized in Table 2. 

The analysis indicates that an increase in chamber pressure from 125 to 200 psia 
will worsen stability only if thrust increased (see Table 2). Improved stability 
is predicted if pc is increased to 200 psia while thrust and injection velocity 
were held constant, irrespective of whether contraction ratio or chamber diameter 
was increased. No appreciable change in stability was predicted for an increase 
in Pc to 200 psia, while holding thrust and fractional injection pressure drop, 
Ap/pc, constant when either contraction ratio or chamber diameter was increased. 

An increase in contraction ratio from 2.0:1 to 3.0:1 improved the predicted sta- 
bility if Pc was allowed to increase so as to maintain constant thrust (Table 2 ) . 
A substantial improvement in predicted stability occurred if the chamber diameter 
was allowed to increase while maintaining constant Pc and thrust. 

A decrease in injection Ap (or injection velocity], while maintaining constant 
thrust, improved the predicted stability when either the number of elements or 
the orifice size was increased. This result is exactly opposite the conclusion 
that might be reached from the Ap curves directly with the original constraints 
on the remaining parameters. The stability decrease caused by the decreased in- 
jection Ap indicated by the curves is more than offset by the stability increase 
caused by the resultant increase in dropsize. 

These results clearly show the importance of carefully defining the manner in 
which parameters are varied, viz., the parameter constraints. Different results 
may be obtained from variations in the same parameter depending on how the remain- 
ing parameters are controlled. Nonetheless, the stability index plots developed 
during this study allow appropriate evaluations to be made. Use of the local 
stability index plots (Fig. B-1 through B-16) is recommended over the averaged 
stability index plots (Fig. B-17 through B-32); the effect of the various param- 
eters on stability appears to be overly suppressed in the latter plots by the 
averaging. Dropsize, and therefore the method used to determine dropsize, has a 
strong effect on predicted stability. Unfortunately, fully satisfactory methods 
of determining dropsize are not available. For this reason the stability index 
plots have been developed for a wide range of dropsizes to that they are not re- 
stricted to particular methods of predicting dropsize. 
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QUASI-LINEAR STABILITY ANALYSES 


The analyses that form the basis for the Northern Research and Engineering Corpo- 
ration (NREC) stability model and, also, the generalized acoustic cavity damping 
model, as well as the previously developed cavity damping model, may be described 
as quasi-linear (i.e., linearized equations are used and nonlinearities are intro- 
duced only through the boundary conditions) . Such analyses are not likely to 
accurately predict strongly nonlinear phenomena but they may be used to reasonably 
represent three-dimensional stability problems. On the other hand, fully nonlinear 
analysis methods are intractable for three-dimensional stability problems. 

The conservation equations for a two-phase mixture may be linearized to obtain an 
inhomogeneous wave equation, as described by Culick (Ref. '42) . The resultant wave 
equation is: 


V 


where 

— R [ 

-r 

This equation was 
small compared to 
is small compared 

bustion. This equation must be solved subject to a boundary condition that may be 
written as: w 

N • Vf * -|p^ It pQ ‘ ^ ‘ ^ 

Culick’ s formulation of several stability problems is based on a first-order 
approximate solution to these equations (Ref. 42). The NREC model is based on the 
same first-order approximate solution with an approximate expression for the com- 
bustion source terms (Ref. 1). The Rocketdyne generalized cavity damping model is 
based on a more accurate iterative-variational solution to these equations with a 
similar expression for the combustion source terms. ' 

The indicated first-order approximate solution may be obtained in the following 
way. A solution to the time- independent homogeneous wave equation and boundary 
conditions corresponding to Eq. 16 and 17 can usually be obtained, i.e.: 

* V *N = 0 


-7 “ -P„ ’ • -7 I? + 


(F - a) 


3t 


. JL 3^ V. - -i- ^ 

2 3t 2 3t 

c c 


(16) 


w + w + Q + 
s s s ^ 




derived under the assumptions that the oscillatory quantities are 
their mean value and, further, that the steady-flow Mach number 
to unity (Ref. 42). The term^represents the oscillatory com- 


N • = 0 


(19) 
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where the subscript N stands for three indices, one for each dimension. Alternate 
boundary conditions can be used, and are desirable in some cases. Assuming a har- 
monic time dependence 3 = u+jot) the time independent forms of Eq.l6 and 17 

may be written as: 

p + k^p = h (20) 


N • Vp = -f 


( 21 ) 


where h and f represent the time independent form of the right-hand sides of Eq. 
16 and 17. If Eq. 18 and 20 are multipled by p and c|)]sj, respectively, then one 
equation is subtracted from the other and the result integrated over the chamber 
volume, the result may be expressed as 


k 


2 




( 22 ) 


The integrals may be evaluated with <()f^ being used in the expressions for h and f 
as an approximation for the pressure 


The latter equation may also be obtained through use of Green's functions with an 
eigenfunction expansion being used for the Green's function. A more useful ex- 
pression has been obtained by Culick from further evaluation of the integrals in 
Eq. 22, i.e.; 


(k2 - /<J.jj^dv = [in 

r_ 2 f - \ - 

i / w <6,, dV + / w Re 1 u - u f • Vd)„dV 
p_c 7 s J s 1 si 


(23) 


NREC INSTABILITY MODEL* 

The NREC model is based on a quasi- linear solution to the conservation equations 
describing the combustion/flow field. NREC developed this model from an analysis 
formulated by Culick in the manner outlined above. The NREC model comprises two 
computer programs (HLMHLT and REFINE). HLMHLT calculates approximate oscillatory 

*Rocketdyne has developed a capability for using the NREC model as a result of a 
subcontract to the General Electric Company for work on contract F3361S-71-C-1742, 
Augmentor Combustion Stability Investigation. Rocketdyne's principal role has 
concerned analytical modeling at the driving processes. 
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pressure distributions and the corresponding frequencies as solutions to the 
Helmholtz equation (time independent form of the wave equation) and without con- 
sidering combustion, i.e., it solves for (f)j^ and kf^ from Eq.18 and 19, usually 
with admittance- type boundary conditions without considering combustion. REFINE 
calculates the frequencies and damping coefficients through a first"-order correc- 
tion for combustion effects. This model can be used to predict the stability of 
longitudinal and radial modes as well as the tangential modes described by the 
Priem model. In principle the model can also describe an acoustic cavity or liner 
but the description is much less accurate than the cavity damping model. The re- 
quired input to the HLMHLT model includes: 

• Chamber dimensions 

• Average sound velocity and gas density 

• Nozzle acoustic admittance 

• Steady flow Mach ntunber at nozzle entrance 

The output from this program (approximate frequencies, damping coefficients, and 
pressure and velocity distributions) are used in REFINE. In addition, the input 
to REFINE includes; 

• A series of coefficients representing the oscillatory combustion 

• Steady flow velocity distribution 

These combustion coefficients may be calculated from models for the spray combus- 
tion. The REFINE program predicts the damping coefficients for each mode corres- 
ponding to the operating conditions. 

The combustion energy release rate is approximated in the NREC model by an expon- 
entially decaying function (with axial position). Employing this assumption, the 
time averaged energy release rate becomes: 

i . p J_ ” C24) 

^ T 

while the oscillatory energy release rate is represented by: 
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The oscillatory energy content (E) and the oscillatory time delay (t) are related 
to the pressure and velocity fluctuations through following combustion coefficients: 



In the original model formulation, droplet vaporization was assumed to be the con- 
trolling mechanism for driving an instability, with the vaporization rate being 
assumed to depend only on the relative droplet Reynolds number. Recent calcula- 
tions (Ref, 43) have shown that this assumption can lead to substantial errors 
with high-frequency instabilities. 

Method of Evaluating Combustion Coefficients 

The mathematical formulation used in the NREC program results in a linear stability 
model (i.e., the model formulation and solution are independent of the pressure 
amplitude). From experimental data, the stability characteristics of a rocket 
engine are known to be dependent upon the instability pressure amplitude. There- 
fore, to allow for these nonlinearities to some extent and to analyze rocket engines 
using the NREC program, the combustion coefficients, which appear in the oscillatory 
energy release rate equation, are considered functions of the pressure amplitude. 

According to the Rayleigh principle, maximum oscillatory driving effect is. obtained 
when an unsteady combustion process (energy addition) occurs in phase with, and in 
a region of, oscillating pressure. However, the overall combustion process comprises 
numerous localized, and often superimposed, combustion-related processes distri- 
buted throughout the chamber. ' Any of these processes may drive or supress a 
resonance, depending upon relative location and sensitivity to the local pressure 
variation. 

A measure of the extent to which the combustion process can reinforce an acoustic 
oscillation is the response factor which is defined so that it is compatible with 
the Rayleigh criterion and is related to the local perturbations (Ref. 43). The 
nonlinear response factor for an unspecified parameter W is defined as: 



[Re(p/p)]^d(wt) 


where W is any oscillatory source of mass or energy and: 

W = W - W (29) 
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The time averaged rate, W, is given by: 
_ 

W _ f W d(wt) 

W J W 2tt 

o o 0 


( 30 ) 


For linear analysis of combustor instability, the response factor can be used 
directly to couple the combustion process to the gas dynamics. In the linear case, 
where all oscillations are sinusoidal, the relationship 

W/W (p/p) 


fully specifies the coupling. However, when oscillatory pressure levels become 
sufficiently large, often at a relative low amplitude, nonlinear effects (wave 
distortion) become important. Further, shock wave-type behavior in either the 
tangential or longitudinal modes has been observed in which the nonlinearity of 
the gas dynamic oscillations is clearly important. To analyze the wave distortion 
problem, the oscillatory pressure used in the calculation of the response factor 
has been assumed to be described by the form (see Ref. 44): 


Re' 


p/p 




cos (nwt) 


(31) 


where pj^ is defined in terms of the peak-to-peak pressure amplitude by (n -*•“): 


-max 


- P, 


min 


P 



(32) 


For low amplitudes the corresponding wave shape is sinusoidal while for large ampli- 
tudes the pressure wave shape approaches a steep fronted shock-type form. 

Through the use of the response factor defined by Eq. 28, the combustion coeffic- 
ients for the NREC program can be evaluated as a function of the pressure ampli- 
tude (Ref. 45) . 

Evaluation of Combustion Coefficients for Spray Combustion 

Theories of droplet combustion are available that may be used to evaluate the ex- 
tent of coupling between droplet buming rate and local pressure and velocity 
fluctuations. In general, droplet burning is enhanced by increased turbulence 
levels or by periodic directional variations in velocity, because droplets are 
relatively heavy and resist following gas streamlines. The overall process is 
controlled largely by heat transfer to the droplet and mass transfer away from 
the droplet, which are, in turn, influenced by the dynamics of the droplet rela- 
tive to its gaseous environment. The droplet vaporization rate may be expressed 
as (see Ref. 46) : 


w, 


6p zk Nu„ 
= s g H 

VAP ^ n 2 

d pv 


(33) 
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z 


(34) 


Cpv % P 



NUh 


2.0 + 0.6 Pr 


1/3 




1/2 


(35) 


The equations used to represent the spray combustion are actually for a single 
vaporizing material but the properties and energy release profiles that are used 
pertain to bipropellant combustion. 


The droplet mass continuity equation is: 
'^VAP " " dz ^^s “zd^ 


(36) 


Substituting Eq. 33 into Eq. 36 and integrating yields the droplet density as a 
function of location. From this density distribution, the energy release rate per 
unit volume may be obtained, i.e.: 



(37) 


where the time delay and the energy content are: 

2 

p, Dj c 

T = ^ pv 

' 6 z kg NUjj 


(38) 


E = AH, 


w^(z=0) 


COMB 


w 


(39) 


Equations 37, 38, and 39 may be expanded into time averaged and oscillatory 

parts. The first of these is the same as that discussed previously. The latter 
equations become: 


t + T' 
T 

o 



z 

o 

T 


Nu, 


HO 



1 


E + E 



E 
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rs^ 
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1 . ^ 


U p 



z Nu„ 


(40) 


(41) 
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where and refer to steady-state (rather than time averaged) conditions and: 




(46) 


(f) 

By comparing the preceding equations with the NREC oscillatory time delay and 
energy content equations, the combustion coefficients may be identified as: 




where 

p/p 



(48) 
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The heat blockage term, z shown in the equations above is related to the combus- 
tion gas and liquid vapor properties by Eq. 34. This blockage term depends on 
the oscillatory droplet temperature because the vapor pressure (p ) at the droplet 
surface is related to the droplet temperature. For a single dropYet, the heating 
rate may be described by (see Ref. 46) : 


dT^ 



z NUj^ 
‘^pd 


T -T, 
g d 

e^-1 



Assuming that, at steady- state conditions. 


(49) 



= 0 


(50) 


and 

P/Pv 

etc., for the other variables, the time averaged solution to the preceding equa- 
tions gives: 


and the oscillatory solution gives: 




1 - iA,\ 

l.A V 

c ' 



(53) 


? = d Jin p^/d Jin Tj 

Examination of the real part of the oscillatory blockage term indicates that the 
real part of the response factor corresponding to this term will always be negative 
and that its magnitude will increase as the frequency of oscillation increases. 
Combining the above results with 
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and noting that the response factor for the Nusselt manber is always positive 
(Ref. 43) lead to the conclusion that the response factor for the time delay will 
be negative for low frequencies and positive for high frequencies. 

The foregoing equations have been used to calculate the combustion *coefficents. 

Based on information given in Ref. 47, 6 has been assumed to be equal to 1.5. The 
response factor has been evaluated with the velocity being assumed to be described 
by: 

Ug = ^ ^ Pj" cos (nwt-6) (55) 

n 

where 8 is the phase angle between the velocity and pressure. When 0 is zero, the 
velocity and pressure are in phase, and traveling wave properties are simulated. 

For standing acoustic modes, the velocity and pressure are out of phase, which 
corresponds to 0 = 90 degrees. 

The steady-state time delay (zq) has been determined by curvefitting results from 
the steady-state combustion model calculations. 

NREC MODEL CALCULATIONS 

The NREC model was used to predict the stability characteristics corresponding to 
some of the same engine conditions that were analyzed with the Pxiem-type model. 

The combustion coefficients were calculated as described above. No acoustic cav- 
ity effects were considered at this stage. The nozzle admittances were calculated 
separately and input to the program. 

Nozzle admittances were calculated with an adaptation of the computer program de- 
scribed by Bell and Zinn (Ref. 48) which is based on the Crocco theory. Typical 
calculated admittance curves are shown in Fig. 7, which pertain to the first tan- 
gential mode and first tangential/ longitudinally coupled modes. The real part of 
the predicted nozzle admittance is positive at high frequencies but negative at 
low frequencies. The negative real components of the admittance implies a driving 
effect which, presumably arises from a conversion of steady-flow energy into oscil- 
latory energy. The nozzle also influences stability through the steady flow, which 
causes a convection of acoustic energy out of the chamber. Because of this con- 
vective effect, the net effect of the nozzle is expected to be a stabilizing one 
even at low frequencies. 

The importance of allowing for oscillations in droplet temperature was investigated 
to determine whether or not this effect should be included in subsequent calculations. 
The effect arises through the heat blockage factor, z, discussed earlier. Stabil- 
ity calculations were made for two chamber sizes. The results are shown in Fig. 8 , 
which show the variation of predicted logarithmic decrements (damping coefficient 
divided by frequency) with the instability amplitude (peak to peak) . Negative 
logarithmic decrements correspond to growing oscillations or unstable conditions. 
Apparently the importance of this droplet temperature effect is greatest at high 
frequencies (i.e., the 2.0-inch-radius chamber) and at this condition the effort 
is significant. Therefore, the effect was included in all subsequent calculations. 
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ACOUSTIC ADMITTANCE OF NOZZLE, A„ = pc lT/yp 



NONDIMENSIONAL FREQUENCY, uR/c 


Figure 7. Predicted Nozzle Admittance for First Tangential and 
Longitudinally Coupled Modes 
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LOGARITHMIC DECRIMENT, Zira/u 



Figure 8. Predicted Oscillatory Decay Rates Showing the Effect of Droplet Temperature 
Oscillation 





Results from a set o£ stability calculations for the first tangential mode in 
three chamber sizes and with two droplet sizes are shown in Fig. 9 and 10. In 
each case a critical amplitude is shown above which an instability is predicted 
the critical amplitude being that corresponding to a logarithmic decrement of zero. 
This critical amplitude, on a zero-to-peak rather than a peak-to-peak basis, 
corresponds to the Priem stability index. 

Results from the NREC model analysis are summarized in Fig. 11 through 13, which 
show the variation of this critical pressure amplitude (zero-to-peak) with drop- 
size, dropsize ratio, and chamber radius. Corresponding curves from the Priem- 
type analysis are also shown. Some of the results from the NREC model analysis 
show a maximum stability index with increasing dropsize. This maximum arises 
because the time average combustion effeciency drops with increasing dropsize, 
which increases the available energy and worsens stability. The reduction in 
efficiency causes E, the energy content (Eq. 34), to increase, thus making more 
energy available for driving an instability. This effect will not occur if the 
combustion efficiency is 100 percent. It is not considered in the Priem-type 
analysis as usually done. The results show a worsening of stability with increas- 
ing chamber radius and, generally, an improvement with increasing dropsize. 

The comparison of results from the NREC model analysis with those from the Priem- 
type analysis (Fig. 11 through 13) shows substantial differences in the magnitude 
of the critical pressure amplitude corresponding to the stability limit. More- 
over, as noted above, the Priem-type results do not show the maximum that appears 
in the NREC-type results. However, there are also major differences in the 
methods of analysis used in these models and, therefore, substantial differences 
in the results are not surprising. Nonetheless, the results show qualitatively 
similar trends, which is encouraging. 
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LOGARITHMIC DECREMENT, 2na/iii 



FRACTIONAL PRESSURE AMPLITUDE, PEAK TO PEAK, 2p/p 

Figure 9. Predicted Oscillatory Decay Rates for = 100 Microns 
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FRACTIONAL PRESSURE AMPLITUDE, PEAK TO PEAK, 2p/p 


Figure 10. Predicted Oscillatory Decay Rates for Dj, = 400 Microns 
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GENERALIZED CAVITY DAMPING MODEL 


A generalized cavity damping model, based on the quasi-linear approach described 
earlier, has been formulated and the preliminary development completed. With 
this model, approximate solutions are obtained to restrictive forms of the inhomo- 
geneous wave equation and boundary conditions (Eq. 16 and 17). The combustion 
source terms were represented in the manner described previously for the NREC 
model. For the vaporization limited process, the source may be written as: 



where 

^ . This expression replaces all of the source terms shown 

in Eq. 16, ! . e. . ^ ^ = R^/C^. 

Development of this cavity damping model has been approached on an incremental 
basis with individual effects being added one at a time. At present, three kinds 
of effects have been introduced: (1) a uniform steady flow, (2) a nozzle admit- 

tance, and (3) pressure-coupled combustion response. However, the steady flow 
contribution at the acoustic cavity has been neglected. (This is probably small 
for a cavity located adjacent to the injector.) The effects of gradients in the 
steady flow, velocity coupled combustion, droplet drag, and steady flow at the 
cavity must be added at a later time. Therefore, at present the generalized 
cavity model is based on solution of the following inhomogeneous wave equation: 


(V^ 


k^) 


p = 2jkM 


« f - 


ijL 

3z^ 




f p (r, e, z’) 



z-z ' 

X dz' 
~T 


= - S(r; p) 


(57) 


i 6t 

where a time dependence of the form e*^*^ has been assumed (g = w + ja). The first 
two terms on the right-hand side of Eq. 57 correspond to the uniform steady flow in 
the axial direction, where M = u/~. The last expression arises from the pressure- 
coupled combustion source terms, where: 


and 


h = _EL_E ^ 
o - s 

X = u T 
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( 58 ) 


Equation 57 has been solved subject to the boundary conditions 
N • Vp= 0 


at rigid walls, including the injector wall and 


N • Vp” _ 


-jky 


(59) 


at the acoustic cavity interface where y is the specific acoustic admittance of 
the cavity. Thus, the effect of the steady flow at the cavity has been neglected. 
The real part of the specific admittance for the boundary condition at the nozzle 
entrance includes the steady-flow effects and may be written as: 


N • 


-jk 


A ^ 

A p 


L 




(60) 


where A = y P N ♦ u/p c , the specific nozzle admittance. The inhomogeneous 
Helmholtz equation has been solved by conversion to an integral equation which 
was solved approximately by an iterative variational technique. The method of 
solution is similar to that used previously for the no-flow, no-combustion case. 


INTEGRAL EQUATION FORMULATION 


The integral equation is developed in terms of the following Helmholtz Green's 
function G (rfr ) satisfying: 


(V^ +Jc2) G C?|r^) = -6 (?-?q) 

N • VG (r\r^) = 0 

An eigenfunction expansion was used for the Green's function, G. 
of a closed cylinder are employed for this expansion, i.e.: 


. 0 


♦n " ''m'V ^ 


COS m 6 
sin m 6 


(61) 

The normal modes 
(62) 


where N represents the three indices m. Si, q. An expression for the Green’s 
function may be readily obtained in terras of these functions, it being; 
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( 63 ) 




''n (''n^ - k^) 


■n 


A®. #> . 

m mx q 


= Jv 


(r)dV 




TT m?fc 0 

2n m=0 (cos m4) only) 





dr J 


m 
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J (a fl) 


mi 


/ dz cos^ ^ = 


L/2 

L 


q # 0 
q = 0 


Through use of this Green's function, the inhomogeneous wave equation and boundary 
condition may be converted to an integral equation (Ref. 49, page 321). 

p(r) = /cCrlr^j) S(r^; pO dV^ + 

* /G(rlr^) N • V?^(r^) dS^ (64) 

where SCr^; p) represents the inhomogeneous term in the wave equation, the 
N • Vp term is given by the boundary conditions, and the vector sign (-*■) has been 
suppressed on r and r for simplicity. This homogeneous integral is solved by an 
iterative-variational technique, 

ITERATIVE-VARIATIONAL SOLUTION TECHNIQUE 

The solution technique is similar to that used previously (Ref. 50). Frequently 
with equations of this form, the integrals may be evaluated with an initial approx- 
imation for pressure to achieve an improved approximation for the pressure. There- 
fore, an iteration form of the expression is written as; 

~(i+l)(Jj ^ j*G(r|r^) S(r^; p*-^^) dV^ + 

/G(r|r^) N • dS^ (65) 
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Because an eigenfunction expansion Vf as used for the Green's function, evaluation 
of the integrals leads to a series expression for the pressure in terms of the 
same eigenfunctions. Thus, from an approximation for p^ ' , an expression is 

obtained of the form: 

^ cos in <t> (66) 

Jl>q 

Evaluation of the integral expression gives an expression for in terms of 

. The iteration is usually started with the eigenfunction corresponding to 
the mode of interest, i.e.: 

i-,q = A,q 

A.q 9* r.q (67) 


= 

Zq 


1.0 

0.0 


The angular dependence drops out because the boundary conditions are uniform in 
the circumferential direction. At times, it is more convenient (if not necessary) 
to use the alV from a simpler case as the starting point for a more complicated 
case, e.g., by tracking the solution as a parameter is incrementally changed. 

Although it is possible, in principle, to use the iteration equation alone to 
solve for complex angular frequency, and determine stability, it is more efficient 
to develop a separate characteristic or eigenvalue equation. The variational 
approach has been used for this purpose. 

Eigenvalue Equation 

Successful use of a variational method depends to a large extent on the form of 
S(r;^) from the wave equation. The variational method described by Morse and 
Ingard (Ref. 49, page 561) has been'used previously,. By separating a particular 
term from the integral expression for pressure, the following equation is obtained 


P(?) 






f-jkj 
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where 


(68) 



Because of the homogeneous form of the equation, an amplitude coefficient can be 
chosen for p such that the term in parentheses ( ) equals unity. Without changing 
the notation for pressure, this choice of amplitude gives: 

P(r) = - jk / Gj^Crjr^) Y(r^) pCr^DdS^ - /Gj^(r|r^) S(r^;p)dV^ 

(69) 


This expression for ^ clearly approaches the normal mode as y and S (r ; p) 
approach zero. Following Morse and Ingard by analogy, a functional is developed 
by multiplying Eq. 69 by -jk / y(r)p(r)dS and then by / S(r;p) dV and adding the 
resultant expressions to Eq. 70 to obtain: 


= 2 yt^o^ pM^S^ S(r^;^dV^ - 

/pCr^) [-jk p(r^)]dS^ +/[ -jk y(r^) p(r^)l ^Gj^Cr | r^) [-jk y(r^) ; (r^)] dV^ dV + 


r-jk y(r^) p(r^)]/GCrlr^) S(r^;^dV^dS - J S(r;|0 p(r)dV + 


fs(r;p) f Gj^(rlr^)[jk yCr^) p(r^) ] dS^dV + /s(r;p)/G(r]r^) S(r^;p0dV^dV 

By varying this equation, a characteristic equation may be obtained: 

^tAwCriN^ - k^)j = -2 f 6 S(t;J() (h^(r) -jk/cj^Cr |r^) y(r^) p(r^)dS^ - 

S(r^?0dV^ -/s(r;^ 6?-dV + p(r) 6S(r^;p)dV^ 


(71) 


(72) 


If the last two terms cancel, Eq, 72 is a variational implying E^. 65 . However, 
the last t^o terms cancel only for the special cases S = p(r) h e’^'^ and 
-M^ ^ p/3z from the wave equation (Eq. 57) the latter case demanding restrictions 
on 5^(r), The latter term is also usually dropped by order of magnitude arguments. 

If only the term S (r;p) = ^ h e is retained, an equation similar to that 

used previously is obtained (Ref. 50), With an approximation for pressure of the 
form ^ = A^(^J, the variation may be performed to obtain: 

/[-jk y(r) p''^^ (r) ] |p^^^ (r) -/G(r|r^) [-jk y(r^) p^^^(r^)]dS^ 

/G(rlr^) S (r^;p^^^)dV^|dS +/s (r;p^^^)( p*^^^ (r) - 

/G(r|r^) [-jk y(r^) p^^^ (r^) ]dS^ - /c(r j r^) S(r;p^^^ )dV^|dV = 0 
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(73) 



This equation may be rewritten as: 

/C-Jlcyp(«) * Is 

The first terra in this equation is the same as that developed previously to assess 
cavity effects alone (Ref. Equation 74 has the interesting form of a weighted 
average of the residual, - ^(i+1) , at the i^h iteration. The weighting factor 

has surface contribution proportional to the wall admittance and a volume contribu- 
tion proportional to the combustion coefficient. Unfortunately, the restriction on 
S(r;p) is undesirable. 

An alternate derivation of an eigenvalue equation corresponding to Eq. 74 was found 
following Morse and Feshbach (Ref. 51). The dervi at ion, involved the use of adjoint 
solutions which could be specified for S(r;^ = TT e , but not in the more gen- 
eral cases. The additional complexity of adjoint functions and adjoint operators 
led to the abandonment of this approach in favor of a "least-squares” method. 

"LEAST- SQUARES" DERIVATION OF EIGENVALUE EQUATION 

Least* squares variational methods are described by M, Becker in his monograph 
(Ref. SI). The related approach used here differs in that the variational is not 
taken to be positive definite, which is the reason for the quotes on "least-squares." 
Consider the functional: 

I * / w(r) |yp - M p|^ dV (75) 

where w(r) is a weighting function ^and the equation 

Vip = Mp 

= /M(r,r^) p(r^)dV^ 

is a shorthand notation for the integral equation (Eq. 64) . The quantity u is an 
artificial eigenvalue. For arbitrary k, there will be discrete eigenvalues, p(k). 

The desired eigenvalues correspond to inverting u(k) = 1 to obtain k = k(iJ = 1). 

The variation of I is taken to include variations of the eigenvalue p: 

61 = 2 6yjw(r)' p(r); |yp(r) - M^dV + 

2/[yp-M^ {u6p-M6p|dV (76) 

This shows that 6l = 0 when the integral equation (Eq. 64) is satisfied. With a 
trial function the variation of I with respect to y implies 



which a new eigenvalue equation for the cavity problem when y = 1 and w(r) is chosen. 
Rather than use the variation of I with respect to 6?' to determine the coefficients 
^£m iterative equation (Eq. 65) was used. Equation 77 holds for 


- p^^’"^^}dV = 0 

(74) 
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arbitrary M or S(r;p) and for arbitrary weighting functions, w(r). Thus, the 
eigenvalue equation is not unique. Different choices of weighting function will 
yield different eigenvalue estimates and, thus, affect the convergence rate in 
the iterative approach. 


One choice of w(r) that is suggested by Eq.73 is: 

w(r) = h e + t -jky(r®)] 6(r - r®)dV 


(78) 


With this choice, Eq. 77 becomes: 

J[ -jky(r) p^^^ (r) 1 I (r) - (r)| dV + 

JSj^(r;p^^^) jp^^^(r) - (r)j dV = 0 


(79) 


where y has been set equal to unity. Equation 79 has the same form as Eq. 74 but 
this derivation shows that it holds for arbitrary S(r,p). 


FORMULATION IN TERMS OF a 


£q 


(i) 


The equations can be placed in a more convenient form for nximerical solution by 
the series coefficients a of Eq. 66. In terms of these coefficients, 

Eq. 65 becomes: 


where 
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and 
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where denotes the slot width and: 
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Note that R , (X s «■) = Kq„,. The terms in Eq. 80 canjbe associated with the 
C£rrespondirtg terms in Eq. 65 via the coefficients, q, M, JT, etc. The last two 
(h^L) expressions correspond to the combust ion- convection integral term in 
Eq. 57^. These two terms must be dropped in the computer program if the steady 
flow M s: 0 t£ avoid division by zero in the program. Analytically, these terms 
vanish when M = 0.0. 


This eigenvalue equation (Eq. 79) can be rearranged by direct substitution to 
obtain: 



q' 


(i)T(M* 
Xq Aq,A' 



(i-1) 

Jl' q' 


)■ 


0 


(81) 


where the notation ”M = 0" denotes that the steady flow terms have been deleted 
from T^ , in_Eq. 80. The no£zle admittance term was retained. Equation 81 
holds ev6n when M 5* 0, but when M = 0, it can be rewritten as: 



(82) 


The essence of Eq. 82 is the residual aj^^ - \ which is weighted by a^^^ 

and some additional factors. Equation s3 was u^ed as the basis for the corap3ter 
program although the derivation is not rigorous for the case M 0. However, 
various modifications of Eq. 82 were tried successfully, e.g.: 



( 83 ) 
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and 


4 - ^Hq rHa ^Zq J " 


Equation 82 seemed to converge most rapidly. Attempts at "least-squares” deriva- 
tion with a more general weighting function 


I = fyp - M ?} wCr.r^D jyp - M p\ dV dV^ = 0 (85) 

ir I Jr^ 

wCr.r^j) = ^ (^.(r) w.. (|) . (r^) 


did not yield Eq. 82. Equation 85 is perhaps more closely related to the method 
of weighted residuals (Ref. 52, page 10). Nevertheless, Eq. 80 and 82 have yielded 
the desired convergence in that additional iterations do not change the complex 
pressure distribution. These resulting solutions are considered to be very good 
approximations to solutions of the original partial differential (Eq. 57). Various 
factors, such as the truncation of the series expansion, Eq, 82, will affect the 
accuracy of the result. 
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COMPUTER RESULTS 


A computer program was developed to solve the foregoing equations and was adapted 
from the previously developed cavity damping model. The computer subroutine used 
to calculate to acoustic impedance of the acoustic cavity was used without change. 

At present, the program is only capable of analyzing the first- tangential model, 
which was done to facilitate model development. Sufficient progress has been made 
to allow preliminary definition of the characteristics of the model predictions, 
but additional work is needed to generalize the model, fully characterize the pre- 
dictions, and relate the model results to measured stability results. 

The computer program has been checked to a large extent by comparison of computa- 
tional results with analytical results from analytically solvable limiting cases. 
Also, comparisons have been made with results obtained witli the previously developed 
model suppressing the new contributions in the generalized model. The two models 
are necessarily somewhat different because of the volume contribution included in 
the new model. 

As noted previously, the model was developed and checked on an incremental basis. 

Ihe results are discussed accordingly. 

Acoustic Cavity 

A comparison of results from the previous cavity model and from the generalized 
shows that they are essentially identical, although different series expressions 
and methods of calculation have been used. Some typical results are shown in 
Fig. 14, 

As shown in Fig. 14, the mode splitting obtained with the previous model was also 
predicted with the generalized model. With the acoustic cavity, two modes are 
obtained in place of the normal first-tangential mode. The lower frequency mode 
approaches the normal chamber mode as the cavity becomes very small. The upper 
frequency mode approaches a quarterwave mode in the cavity as the cavity width 
becomes small; nonetheless, the chamber is fully involved in this oscillation for 
cavities wide enough to be useful for stabilization. 

The damping associated with each of these modes has not been calculated beyond 
the crossover point between the two damping coefficient curves because only the 
most weakly damped mode is likely to occur. 

For the lower frequency mode, the oscillatory pressure amplitude tends to increase 
with distance from the cavity. For the high frequency mode, the oscillatory pres- 
sure amplitude tends to decrease with distance from the cavity. It is of interest 
to note these pressure distributions could be approximated by only the first two 
terms in the series expansion, i.e.: 

P =Jj r/r^) cos m9 . 8,5^ cos (86) 
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It is also of interest that the flow of energy can be related to phase of the 
pressure: 


The energy flux is (.for no steady flow) 

I = p Re (u) 

so that the time averaged energy flux is: 


( 87 ) 


( 88 ) 


<T> = <p Re{u}> = - I:- Vtp 

pa) 


(89) 


The lines of constant phase are expected to have the qualitative features sketched 
below: 



An examination of the computer results shows that this is the case if spatial oscil- 
lations in are averaged out, since spatial oscillations are probably due to the 
finite nature of the expansion series. Physically, acoustical energy flows out 
of the chamber at the acoustic cavity. This energy flows from all parts of the 
chamber. 

Uniform Combustion 

The simplest source term is a spatially uniform combustion. For this case the 
wave equation becomes : 

(V^ + k^) 5" = _hp (90) 


For the rigid wall (no cavity case) , 


P = 


100 


jkct 

e 


the solution to this equation is: 


k 




2n 


h 

100 


(91) 
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For a purely imaginary combustion coefficient, a growing wave results; the com- 
puter program was checked for this simple case. 

Spatially Uniform Combustion With Acoustic Cavity 

A set of calculations was made for a uniform combustion source term. Results from 
these calculations are shown in Fig. 15. The primary effect is to shift the damp- 
ing coefficient curve downward from the no-source case, as expected. Also, the 
net damping coefficient approximately equals the sum of the cavity and source 
contributions (calculated individually), i.e.: 

*^comb * "^cav %oth ^ ^ 

An examination of the calculated pressure distributions shows little change due 
to the addition of the uniform combustion. 


Nonuniform Combustion With Acoustic Cavity 

Calculations were made with an exponentially varying combustion distribution, i.e. 

h(z) = h e C93) 

From Eq. 70 it is evident that the contribution to the overall damping coefficient 
should be approximately proportional to: 


Ira Jh(z)} dV 


(94) 


This expression suggests that this contribution will be nearly constant as X is 
varied if p is approximately independent of axial position and: 

Jim |h(z)| dz = constant (95) 


or 


h(X) = h(X - “) X — ^/X 

1 - e 


(96) 


The effectiveness of this approach is demonstrated by the following results (no 
acoustic cavity) ; 


, inches 

h 

o 

X, inches 

cor 

w 

c 

ar 
w 

c 

1. 

45 

0.0 + 0.006j 

1000 

1.8413 

-0.02453 



0.0 + 0.013878j 

5 

1.8414 

-0.02466 



0.0 + 0.0018783J 

-5 

1.8414 

-0.02466 
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DAMPING COEFFICIENT, a, SEC 





The change in damping (or driving) is indicated by the change in the The 

negative value for X corresponds to shifting the cavity to the opposite”end of 
the chamber. In this mode, there is a small average flow of acoustic energy in 
the 2-direction from the region of greater combustion to that of less combustion. 
When an acoustic cavity was introduced, the following results were obtained: 



s 

Re 


Imjh^f 

X, inches 

o)r 

w '■ 

c 

ar 
w 

c 

1.45 

0, 

.0 

0.013878 

5.0 

1.7756 

0.02116 





0.006000 

1000.0 

1.7764 

0.02331 





0.001878 

-5.0 

1.7767 

0.02488 

1. 

65 



0.013878 

5.0 

1.8819 

0.02066 





0.006000 

1000.0 

1.8832 

0.01901 





0.001878 

-5.0 

1.8845 

0.01732 


Thus, the effect of changing X has an opposite effect for the lower and upper 
(£5 = 1.45, 1.65 correspond to the lower and upper modes, respectively). This 
result implies that the effects are not additive, i.e.: 

ct (h, y ,) ?^ a (h = 0, y„„„) + a (h, y =0) 

'' ’ ''cav^ ^ * •'cav^ '• * •'cav ^ 

The trends of these results can be explained in terms of the following equation: 

a ~ J Imjh(z)} IpI^ dV (97) 

The damping due to the cavity is greatest when the trends in pressure amplitude 
and combustion source, h(z), are the same, i.e., both increasing or both decreas- 
ing with axial position. 

Effects of Steady Flow 

At present the model only includes a uniform steady flow, M = u/c, independent of 
position. Steady flow affects stability through its convective effect in the wave 
equation (Eq. 57) and also through the nozzle boundary condition (Eq. 60). Also, 
the injector boundary condition has been maintained as 3^/3z= 0. 

Without combustion or acoustic cavities, the inhomogeneous wave equation can be 
solved analytically. However, the first-tangential mode is independent of z and 
thus the addition of steady flow has no effect, since the relevant terms involve 
z-derivatives of pressure. Therefore, the program was checked out by comparison 
with analytical results for the first-tangential/first- longitudinal mode <J»201 " 
Jl(1.841 r/R) cos ttz/L cos 9, with injector and nozzle boundary conditions of 
99'/3n = 0. The results were in agreement. 

Calculations to investigate the effects of steady flow with distributed combustion 
but no cavity, with 

hp = 0.0 +0.013878j 
X =5.0 inches 
M = 0.1 
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led to the conclusion that there was no significant difference between the flow 
and no-flow cases. However, when an acoustic cavity was also included, the effects 
were significant, as shown below: 


X, inches 


i , inches 
s* 

h 

o 

1.45 

0.0 + 0.01387824j 
0.0 + 0.01387824j 
0.0 + O.Oj 
0.0 + O.Oj 



(jor 
w 

ar 

w 

M 

c 

c 

1.7751 

0.01729 

0.0 

1.7798 

0.04511 

0.1 

1.7707 

0.04280 

0.0 

1.7753 

0.07190 

0.1 


The change in the damping coefficient, a, due to a steady-state flow M = 0.1 
is approximately independent of the presence of combustion r^f/c = 

0.045 - 0.017 = 0.028 and 0.072 - 0.043 = 0.029). The presence of steady flow 
increased the damping rate a for the cases considered, but for these cases, the 
amplitude of the oscillatory pressurelpTl decreases with increasing axial position. 

Effects of Nozzle Admittance 


For these calculations, the nozzle admittance has been taken to be a purely real ■ 
positive number in Eq. 82 and, thus, has a dan 5 >ing effect. The major effect of 
nozzle admittance is to chainge the phase of the pressure distribution so as to 
imply flow of acoustic energy into the nozzle according to Eq. 97. This can be 
seen in the following case (no acoustic cavity, X = 1000 inches) : 


X, inches 


0.0 + 0.006j 
0.0 + 0.006j 


1.8413 

1.8432 


•0.02453 

-0.00849 


The phase angle can be shown to be a function only of z in the absence of the 
acoustic cavity and is shown in Fig. .' 16. These results agreed with analytical 
results for the same case. With an acoustic cavity, the steady-flow effect was 
changed: 





car 

ar 

X, inches 

h 

i, , inches 

w 

w 


0 

s 

c 

c 

5 

0.0 + 0.01387824j 

1.45 

1.7798 

0.04511 

5 

0.0 + 0.01387824j 

1.45 

1.7735 

0.04922 


The nozzle has the expected damping effect but, although the admittance is the 
same as in the preceding case, the incremental change in damping (Aar^/c = 0.04922 
to 0.04511 = 0.00411) is down by a factor of four from the preceding case fAar^^/c = 
-0.00849 + 0.02453 = 0.01604). It is noted that the pressure amplitude Ip at the 
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nozzle has the same radial distribution in both cases, but the pressure |p[ decreases 
by approximately 50 percent in going from the injector to the nozzle in the latter 
case. The factor of four is most likely determined by the square pressure of the 
pressure at the nozzle in these cases. 

It is expected that nozzle admittance will more strongly damp the lower mode because 
the pressure at the nozzle is relatively higher. The plot of phase angle in Fig. 16 
shows that the flow of energy to the nozzle in accord with Eq. 97. The phase depends 
strongly on r for 0 < z/L <0.2 because of the acoustic cavity and, therefore, is 
not shown. 

Effect of Cumulative Combustion-Convection Term 

The last source term in Eq. 57 may be described as the cumulative combustion- 
convection term. The inclusion of th^s term in the program was found to have 
negligible effect on the results for M = 0.1 for the few cases studied. The 
reason is probably associated with the rapid oscillations corresponding to the 
exponential term in the integral. 
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CONCLUDING REMARKS 

The principal objective of the program described herein was to analytically pre- 
dict the relative combustion stability of various propellant combinations of 
interest to NASA-JSC when used with hardware configurations representative of 
current design and properties with or without acoustic cavities. As noted earlier, 
consideration was Originally given to propellant combinations of the LOX/hydrocarbon, 
LOX/amine, and NTO/araine families, but was later concentrated on the NTO/amine fam- 
ily because of diminished interest in the other families, 

To meet the objective of the program, several kinds of analysis were performed. Al- 
though this objective was achieved within the context of the contracted program and 
considerable progress has been made, the results show the need for further analysis. 

Considerable information was generated as a result of the literature review. How- 
ever, the information shows considerable diversity, which probably reflects the gen- 
eral complexity of the instability problem. The results are probably best utilized 
in conjunction with an analytical model to relate the conditions at a point of in- 
terest to those that previously evaluated experimentally. 

A Priem-type stability analysis was done for a relatively broad range of conditions. 
Extensive plots of the variation of the stability index, Ap, with important physical 
parameters have been developed. Because these plots were developed in terms of 
physical parameters, which minimized the total number of parameters, they do not 
show all effects of interest. However, the plots may be used by inteirpolation to 
assess the effect of varying other parameters and, similarly, they may be used to 
assess the stability of new hardware configurations with minor additional calcula- 
tions. The manner in which this may be done has been illustrated by the additional 
calculations. The results from this analysis show the importance, when performing 
stability analyses, of the choice of parameters to be held constant. Opposite 
effects can be predicted by changing the parameters held constant. In addition, 
the results suggest that the propellant mass flux within an engine has a very 
strong influence on stability. The greatest stability changes were predicted when 
Cl) chamber pressure was increased by increasing the propellant flow through a 
fixed set of hardware (worsened stability) or (2) contraction ratio was increased 
with a fixed thrust (improved stability) . 

The analysis done with the NREC stability model was less extensive than that done 
with the Priem-type model. The Priem-type analysis was done with generalized pa- 
rametric curves that were available; but similar curves are not available for the 
NREC analysis and greater effort was required for each of these cases. The re- 
sults from the NREC analysis appear to be roughly compatible with those from the 
Priem-type analysis. The degree of agreement appears consistent with the major 
differences in approach that have been used in development of the two models. The 
NREC model, or similar approaches, appears to be a valuable method of analysis for 
liquid-propellant engines. The Priem and NREC models complement each other because 
each includes important factors not included in the other. Further work with the 
NREC model is recommended. 
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Because o£ the need for a model to aid in the design of acoustic cavities that 
specifically included the effects of combustion and steady flow, development of 
a generalized acoustic cavity model was also undertaken, This model was based on 
a combination of the concepts used in the existing cavity damping model (which 
does not specifically allow for combustion and steady flow) and those used in the 
NREC model. Preliminary development of this model was completed, with an iterative 
variational solution method being used. At present the model includes effects due 
to pressure- coup led combustion response (velocity coupling effects remain to be 
added), nozzle effects, and steady flow (as a uniform approximattion) , as well as 
the acoustic cavity. Each of these contributions has been added incrementally and 
checked by comparison with analytical results for limiting cases. The limited re- 
sults obtained thus far show the importance of interactions between the effects of 
the cavities and those due to the nozzle, combustion and steady flow, under some 
circumstances. Because of these interactions, the need to adequately allow for 
them in cavity design is evident and additional work is recommended. 


R-9353 

62 



REFERENCES 


1. Dix, M. D. and G. E. Smith, "Analysis of Combustion Instability in Aircraft 
Engine Augmentors," AIAA Paper No. 71-700, presented at the AIAA/SAE 7th 
Propulsion Joint Specialist Conference, Salt Lake City, Utah, June 1971. 

2. Culick, F. E. C., "Stability of High-Frequency Pressure Oscillations in 
Rocket Combustion Chambers," AIAA J ., Vol . 1, No. 5, May 1963. 

3. Bates, J. 0., Experimental Evaluation of Storable Liquid Propellants in 
Large Rocket Engines , R-1955, Rocketdyne Division, Rockwell International, 
Canoga Park, California, December 1959. 

4. Levine, R. S., Rocketdyne Experience Regarding the Combustion Chamber Tract- 
ability of Various Propellant Combinations, RR-58-45, Rocketdyne, November 
1958. 

5. Nestlerode, J., Rocketdyne, Private Communication. 

6. Mikuni, D., Rocketdyne, Private Communication. 

7. Stassinos, J., Rocketdyne, Private Communication. 

8. Rocketdyne Research Memorandum, RM 848/352, 1 February 1962, 

9. Coultas, T. A. and R. S. Levine, "Baffles and Chemical Additives as Acoustic 
Combustion Stability Suppressors in a Two-Dimensional Thrust Chamber," Second 
ICRPG Combustion Conference, CPIA Publication No. 105, Vol. 2, May 1966, p. 61. 

10. Rocketdyne Research Memorandum, RM 413/91, 18 June 1959.' 

11. Fourteenth NASA F-1 Program Review, Rocketdyne, May 1964. 

12. Rocketdyne Research Memorandum, CHTUM 58-31, 25 August 1958. 

13. Richtenburg, W,, Rocketdyne, Private Communication, 

14. Shuster, E., Rocketdyne, Private Communication. 

15. Oberg, C. L., "LM Ascent Engine Acoustic Cavity Study," Sixth ICRPG Combus- 
tion Conference, CPIA Publ. No, 192, Vol. 1, December 1969, p. 303. 

16. Oberg, C, L. et al., "Evaluation of Acoustic Cavities for Combustion Stabili- 
zation," Seventh JANNAF Combustion Meeting, CPIA Publ. No. 204, Vol. 1, 

February 1971, p. 743. 

17. Oberg, C. L, et al . , "Combustion Stabilization with Acoustic Absorbers," 

Fifth ICRPG Combustion Conference, CPIA Publ. No, 183, December 1968, p. 359. 

18. Bell Aerosystems Report No. D8539-953001 . 


R-9353 

63 



19. Senneff, J. M. and P. J. Morgante, "Combustion Stability Investigation of 

the LM Ascent Engine," Second ICRPG Combustion Conference, CPIA Publ. No. lOS, 
Vol. 1, May 1966, p. 23. 

20. Senneff, J. and K. Berman, "Stability Investigation Relating to Bell/LM Ascent 
Engine," Sixth ICRPG Combustion Conference, CPIA Publ. No. 192, Vol. 1, 
December 1969, p. 317. 

21. Sherman, E. W., Jr. and D. T. Harrje, "History of Model 8533 Evolution with 
Particular Emphasis on Stability Aspects,: Sixth ICRPG Combustion Conference, 
CPIA Publ. No. 192, Vol. 1, December 1969, p. 351. 

22. Reardon, F. H., "Advanced Agena Engine Combustion Stability Program,: Sixth 
ICRPG Combustion Conference, CPIA Publ. No. 192, Vol. 1, December 1969, p. 357, 

23. Mehegan, P., Rocketdyne, Private Communication. 

24. Dykema Owen, Aerospace, Private Communication. 

25. Hefner, R. J., "Combustion Stability Development at Aerojet-General, First 
ICRPG Combustion Instability Conference, CPIA Publ. No. 68, Vol. 1, January 
1965, p. 9. 

26. Aerojet-General Corporation Report GEMSIP FR-1, SSD-TR-66-2, 31 August 1965. 

27. McBride, J. M. and R. J. Hefner, "Combustion Stability Characteristics of the 
Apollo SPS and Transtage Liquid Rocket Engines," Fourth ICRPG Combustion Con- 
ference, CPIA Publ. No. 162, Vol. 1, December 1967, p. 29. 

28. Abbe, C. J. and W. L. Putz, "The Bipropellant Combustion Stability Character- 
istics of Hydrazine and MMH Fuels with Nitrogen Tetroxide for Post-Boost Pro- 
pulsion Application," Sixth ICRPG Combustion Conference, CPIA Publ. No. 192, 
Vol. 1, December 1969, p. 93. 

29. Abbe, C. J. et al . , "Influence of Storable Propellant Liquid Rocket Design 
Parameters on Combustion Instability," Fourth ICRPG Combustion Conference, 

CPIA Publ. No. 162, Vol. 1, December 1967, p. 73. Also see J. Spacecraft and 
Rockets, Vol. 5, No. 5, May 1968, pp. 584-90. 

30. Weiss, R. and R. Klopotek, Experimental Evaluation of the Titan III Transtage 
Engine Combustion Stability Characteristics , AFRPL Report No. TR-66-51, 

March 1966. 

31. Houseman, J., "Jet Separation and Popping with Hypergolic Propellants," 

Seventh JANNAF Combustion Meeting, CPIA Publ. No. 204, Vol. 1, February 1971, 
p. 445. 

32. Clayton, R. M., "Baffle Performance with High- Impedance Injectors," Seventh 
JANNAF Combustion Meeting, CPIA Publ, No. 204, Vol. 1, February 1971, p. 797. 


R-9353 

64 



33. Baileyj C. R., "Acoustic Liner for the C-1 Engine," Fifth ICRPG Combustion 
Conference, CPIA Publ. No. 183, December 1968, p, 335. 

34. Hoeptner, H. W., "Test Results and Analysis of Pulse Initiated Oscillations 

in the UTC 8K Engine," Second ICRPG Combustion Conference, CPIA Publ. No. 105, 
Vol. 1, May 1966, p. 47. 

35. Moberg, D. A., "Attentuation of Tangential Combustion Instability in an 
Ablative, Hypergolic Bipropellant 1750 -Pounds -Thrust Rocket Engine," Fourth 
ICRPG Combustion Conference, CPIA Publ. No, 162, Vol. 1, December 1967, p. 15. 

36. Minton, S. J. and E. Z. Zwick, "An Investigation of Manifold Explosions in 
Rocket Engines," Third ICRPG Combustion Conference, CPIA Publ. No. 138, Vol. 1, 
February 1967, p. 497. 

37. Marquardt Report S-483, An Investigation of Oxidizer Manifold Explosions in 
the Apollo SM/RCS Engine, 1 February 1966. 

38. Coultas, T. A. and R. C. Kesselring, "Extension of the Priem Theory and Its 
Use in Simulation of Instability on the Computer," Proceedings of the Second 
ICRPG Combustion Conference, CPIA Publication No. 105 Chemical Propulsion 
Information Agency, Silver Spring Maryland, May 1966, page 163-192. 

39. Lambiris, S., L. P. Combs, and R. S. Levine, "Stable Combustion Processes in 
Liquid Propellant Rocket Engines," Combustion and Propulsion, Fifth AGARD 
Colloquium (MacMillan Co., New York, 1963), pages 569-634. 

40. Ingebo, R. E., Dropsjze Distributions for Impinging Jet Breakup in Airstreams 
Simulating the Velocity Conditions in Rocket Combustors , NASA TN-4222, National 
Aeronautics and Space Administration, Washington, D.C., 1958. 

41. Priem, R. J. and D. C. Guentert, Combustion Instability Limits Determined by 
a Non-Linear Theory and a One-Dimensional Model, NASA TN D-1409, National 
Aeronautics and Space Administration, Washington, D.C., October 1962. 

42. Culick, F. E. C., Interactions Between the Flow Field, Combustion and Wave 
Motions in Rocket Motors. NWC TP-5349, Naval Weapons Center, China Lake, 
California, June 1972. 

43. Heidmann, M.F. , Amplification by Wave Distortion of the Dynamic Response of 
Vaporization Limited Combustion , NASA TND-6287, National Aeronautics and 
Space Administration, Washington, D.C., May 1971. 

44. Heidmann, M. F., Empirical Characterization of Some Pressure Wave Shapes in 
Strong Travelling Transverse Acoustic Modes, NASA TMX-1716, National Aero- 
nautics and Space Administration, Washington, D.C., January 1969. 

45. Schuman, M. D. and C. L. Oberg, Rocketdyne contribution to final report 
Augmentor Combustion Instability Investigation, General Electric Company, 
Evandale, Ohio, June 1974. 


R-9353 

65 


46. Priem, R. J. and M. F. Heidmann, Propellant Vaporization as a Design Criterion 
for Rocket Engine Combustion Chambers , NASA Technical Report R-67, National 
Aeronautics and Space Administration, Washington, D.C. 1960. 

47. Heidmann, M. F. and P. R. Wieber, Analysis of Frequency Response Character- 
istics of Propellant Vaporization, NASA TND-3749, December 1966. 

48. Bell, W. A. and B. T. Zinn, The Prediction of Three-Dimensional Liquid- 
Propellant Rocket Nozzle Admittances, NASA CR- 12 1129, Georgia Institute of 
Technology, Atlanta, Georgia, February 1973. 

49. Morse, P. M. and K. U. Ingard, "Theoretical Acoustics," McGraw-Hill , New York 
(1968). 

50. Oberg, C. L. et al.. Final Report, Evaluation of Acoustic Cavities for Com- 
bustion Stabilization, NASA CR-115007, R-8757, Rocketdyne, 1971. 

51. Morse, P. M. and H. Feshback, Methods of Theoretical Physics, McGraw-Hill , 

New York (1953), Vol. I, p. 1108. 

52. Becker, M., The Principles and Applications of Variational Methods , MIT 
Press, Cambridge, Massachusetts (1964). 


R-9353 

66 



APPENDIX A 

DROPSIZE CORRELATIONS 

A number o£ correlations have been developed for predicting propellant drops ize. 
Some of these are summarized below. 

• Ingebo (Ref. A-1) for like doublets: 



where 


K = 0 




1/4 


= 10.7 X lO"^ g/cc 
P = Pj,/RT 


V ~ fps 
d ~ inch 


|AV| ~ V . . D_.-^ microns 

inj 30 

D^q = mass mean drop diameter 

The velocity difference [AvI = |Vg = | may be approximated by neglect 

ing Vg or, when used with a steady-state"^ combustion model, it may be ad- 
justed until the predicted performance equal the measured performance 


Dickerson (Ref. A-2) for like doublets: 

,0.568 


D = 5.84 X 10 


4 d" 
V 


0.852 


(PCFj 


where 


D = 2.2 D^q ~ microns 
d ~ inches 
V ~ ft/sec 
PCF = 1.54 


Ib/ft sec 
^ dyne/ cm 
~g/cc 
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Zajac (Ref. A- 3) for like doublets (60-degree impingment angle)^ 

A n p “0-52 

Laminar Jet D = 4.85 x 10 V ' d 0.57 

wax Pj 

— 4 -1 0 0 57 Pc 

Turbulent Jet D = 15.9 x 10 V d^'^ (^) 

wax rj 

where 


D ^ microns 
V ~ ft/sec 
d ~ inches 

Pc/Pj = Vel. profile parameter ~1 

and 

= (1.42 - .0073y)D^q 

where 


also 


= dropsize for 60=degree impingement angle 
Y = impingement angle 


D = 


D 

wax 


K 


where 


1/4 
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APPENDIX B 


RESULTS FROM PRIEM-TYPE STABILITY ANALYSIS 

This appendix contains a series of plots of the calculated stability index, Ap 
for a variety of conditions. These plots summarize the results from the Priem 
type stability analysis. 
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Fuel Dropsize, D^, microns 

Figure B-1. Local Stability Index as a Function of Dropsize 

for P = 125 psia, AP. . = 0.2 P , CR = 2 
c inj c* 
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Figure B-2 


Fuel Dropsize, D^, microns 

Local Stability Index as a Function o£ Dropsize 

for P = 125 psia, AP. . = 0.2 P , CR = 3 
c inj c 
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Figure B-3. Local Stability Index as a Function of Dropsize 
for P = 125 psia, AP. . = 0.4 P , CR = 2 
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Figure B- 


Fuel Dropsize, D^, microns 

. Local Stability Index as a Function of Dropsize 
for = 125 psia, =0.2 P^, CR = 2 
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Fuel Dropsize, D£, microns 

Figure B-5. Local Stability Index as a Function of Dropsize 

for P = 200 psia, AP, . = 0.2 P , CR = 2 
c * inj c 
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Fuel Dropsize, D^, microns 

Figure B-7. Local Stability Index as a Function of Dropsize 
for = 200 psia, = 0-4 P^., CR = 2 


R-9353 




Stability Index, 



Figure B-8. Local Stability Index as a Function of Dropsize 
for = 200 psia, = 0.4 P^, CR = 3 
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Fuel Dropsize, D^, microns 

Figure B-9. Local Stability Index as a Function of Dropsize 
for = 125 psia, =0.2 P^, CR = 2 
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Figur© B-10. Local Stability Indox as a Function of Drops 
for = 125 psia, AP_ = 0.2 P^, CR = 3 
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Fuel Drops ize, microns 

Figure B-14. Local Stability Index as a Function of Drops! ze 
for P = 200 psia, AP. - = 0.2 P , CR = 3 

^ ITll C. 
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Stability Index, 
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Fuel Dropsize, D^, microns 

Figure B-16. Local Stability Index as a Function of Dropsize 

for P = 200 psia, AP. . * 0.4 P , CR = 3 
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Fuel Dropsize, D^, microns 

Figure B-17. Cumulative Stability Index as a Function of Dropsize 

for = 125 psia, AP. - = 0.2 P^, CR = 2 
c inj c' 
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Figure B-18. Cumulative Stability Index as a Function of Dropsize 
for P = 125 psia, AP. . = 0.2 P , CR = 2 
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Fuel Drop size, D£, microns 

Figure B-19, Cumulative Stability Index as a Function of Dropsize 

for P = 125 psia, AP. . = 0.4 P , CR = 2 
c inj c 
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Fuel Dropsize, microns 

Figure B-20. Cumulative Stability Index as a Function of Drop; 
for Pj. = 125 psia, AP_ = 0.2 P^, CR = 3 
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Fuel Drops ize, D^, microns 

Figure B-21. Cumulative Stability Index as a Function of Dropsize 
for P = 200 psia, AP. . = 0.2 P CR = 2 

C J ^ 
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Fuel Dropsize, D^, microns 

Figure B-22. Cumulative Stability Index as a Function of Drops 
for = 200 psia, = 0.2 P^, CR = 3 
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Fuel Drops ize, D^, microns 

Figure B-23. Cumulative Stability Index as a Function of Dropsize 
for = 200 psia, AP^nj = 0.4 P^, CR = 2 
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Stability Index, 


= 200 psia 

AP. . = 0.4 P 
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Figure B-24. Cumulative Stability Index as a Function of Dropsize 
for P^ = 200 psia, AP_ = 0.4 P^, CR = 3 
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Fuel Dropsize, D^, microns 

Figure B-25. Cumulative Stability Index as a Function of Dropsize 

for P = 125 psia, AP. . = 0.2 P^, CR = 3 
c inj c’ 
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Stability Index, 
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AP. . * 0.2 P 
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Fuel Dropsize, D^, microns 


Figure B-26. Cumulative Stability Index as a Function of Dropsize 
for = 125 psia, = G.2 P^, CR - 2 
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Stability Index, 




Fuel Dropsize, D^, microns 

Figure B-27. Cumulative Stability Index as a Function of Dropsize 
for P = 125 psia, AP. . = 0.4 P , CR = 2 
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Fuel Dropsize, microns 

Figure B-28. Cumulative Stability Index as a Function of Dropsize 

for P = 125 psia, AP. . = 0.4 P , CR = 3 
c inj c 
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Figure E-29. Cumulative Stability Index as a Function of Dropsize 

for P = 200 psia, AP. . = 0.2 P , CR = 2 
c inj c' 


R-9:353 

100 




Stability Index, 


= 200 psia 

= 0.2 P 

nj c 

CR = 3 



Fuel Dropsize, D^, microns 


Figure B-30. Cumulative Stability Index as a Function of Dropsize 

for P = 200, AP. . = 0.2 P , CR = 3 
c in j c 
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Fuel Dropsize, D^, microns 

Figure B-31. Cumulative Stability Index as a Function of Dropsize 

for P = 200 psia, &P. . = 0.4 P , CR = 2 
c inj c ’ 
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Fuel Dropsize, D^, microns 


Figure B-32. Cumulative Stability Index as a Function of Dropsize 

for P = 200 psia, AP. . = 0.4 P , CR = 3 
c in j c 
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